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Abstract 

In the first part of tfiis review we discuss tfie basic observational features at the end of 
the cosmic ray energy spectrum. We also present there the main characteristics of each 
of the experiments involved in the detection of these particles. We then briefly discuss 
the status of the chemical composition and the distribution of arrival directions of 
cosmic rays. After that, we examine the energy losses during propagation, introducing 
the Greisen-Zaptsepin-Kuzmin (GZK) cutoff, and discuss the level of confidence with 
which each experiment have detected particles beyond the GZK energy limit. In the 
second part of the review, we discuss astrophysical environments able to accelerate 
particles up to such high energies, including active galactic nuclei, large scale galactic 
wind termination shocks, relativistic jets and hot-spots of Fanaroff- Riley radiogalaxies, 
pulsars, magnetars, quasar remnants, starbursts, colliding galaxies, and gamma ray 
burst fireballs. In the third part of the review we provide a brief summary of scenarios 
which try to explain the super-GZK events with the help of new physics beyond the 
standard model. In the last section, we give an overview on neutrino telescopes and 
existing limits on the energy spectrum and discuss some of the prospects for a new 
(multi-particle) astronomy. Finally, we outline how extraterrestrial neutrino fluxes can 
be used to probe new physics beyond the electroweak scale. 

Solicited Review Article Prepared for Reports on Progress in Physics 
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PREFACE 



Reviewing cosmic ray physics is a risky business. Theoretical models continue to 
appear at an amazing rate, both in the astrophysical and more exotic domains. We 
warn the reader: we do not (nor we could) intend to make an homogeneous coverage of 
all the ideas of our fellow colleagues. Reading differently focused reviews on the issue 
(quoted here along the way) is, in our view, the best approach to such a wide topic of 
research. 
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I. THERE'S SOMETHING ABOUT COSMIC RAY OBSERVATIONS 



A. Experiments and future projects 



The cosmic ray (CR) spectrum spans over roughly 11 decades of energy. Contin- 
uously running monitoring using sophisticated equipment on high altitude balloons 
and ingenious installations on the Earth's surface encompass a plummeting flux that 
goes down from 10^ m~^ s~^ at ~ 10^ eV to 10~^ km^^ yr^^ at ~ 10^° eV. Its shape 
is remarkably featureless, with little deviation from a constant power law across this 
large energy range. The small change in slope, from oc E~'^''^ to oc E~^'^, near 10^^'^ eV 
is known as the "knee". The spectrum steepens further to E'^'"^ above the "dip" 
(fs lO^^-'^ eV), and then flattens to E'^-^ at the "ankle" 10^^ eV). Within statistical 
uncertainty of current observations, which is large above 10^° eV, the upper end of 
the spectrum is consistent with a simple extrapolation at that slope to the highest 
energies, possibly with a slight accumulation around 10^^'^ eV (For recent surveys of 
experimental data the reader is referred to [H Q, 0, 0, 

It is a lucky coincidence that at the energy (~ 10^^ eV) where direct measurement 
of CRs becomes limited by detector area and exposure time, the resulting air showers 
that such particles produce when they strike the upper atmosphere become big enough 
to be detectable at ground level. There are several techniques that can be employed 
in the process of detection: 

(i) Direct detection of shower particles is the most commonly used method, and 
involves constructing an array of sensors spread over a large area to sample particle 
densities as the shower arrives at the Earth's surface. The pioneering development of 
the air shower techniques (and the first use of plastic scintillation detectors for the dual 
use of measuring arrival directions and particle densities) was started at the Agassiz 
Station of the Harvard College Observatory, a work carried out between 1954 and 
1957 0, 0, 01- The existence of primary particles with energies greater than 10^^ eV 
was established by the observation of one shower with more than 10^ particles. Soon 
afterwards, this technique flourished with measurements of ultra high energy cosmic 
rays (UHECRs) with the Volcano Ranch experiment in the 60 s [13, IllL Il2l | , as well as 



with ^veral other array., =„ch a. Haverah Park in England 0, Snt=k in Ru.=ia Q 

0], the Sydney University Giant Airshower Recorder in Australia (SUGAR) and 
the Akeno Giant Air Shower Array (AGASA) in Japan jl7. 18|. 

(ii) Another well-established method of detection involves measurement of the lon- 
gitudinal development (number of particles versus atmospheric depth) of the extensive 
air shower (EAS) by sensing the fluorescence light produced via interactions of the 
charged particles in the atmosphere. The emitted light is typically in the 300 - 400 nm 
ultraviolet range to which the atmosphere is quite transparent. Under favorable at- 
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mospheric conditions, EASs can be detected at distances as large as 20 km, about 2 
attenuation lengths in a standard desert atmosphere at ground level. However, obser- 
vations can only be done on clear Moonless nights, resulting in an average 10% duty 
cycle. The fluorescence technique has so far been implemented only in the Dugway 
desert (Utah). Following a successful trial at Volcano Ranch [3| the group from the 
University of Utah built a device containing two separated Fly's Eyes (20l |21|. The 



two-eye configuration monitored the sky from 1986 until 1993. As an up-scaled ver- 
sion of Fly's Eve , th e High Resolution (HiRes) Fly's Eye detector begun operations 



in May 1997 [2j,|23|. In monocular mode, the effective acceptance of this instrument 
is ~ 350(1000) km^ sr at 10^^ (10^°) eV, on average about 6 times the Fly's Eye ac- 
ceptance, and the threshold energy is 10^^ eV. This takes into account a 10% duty 
cycle. 

(Hi) A more recently proposed technique uses radar echos from the column of ionized 
air produced by the shower. This idea suggested already in 1940 j^l, has been recently 
re-explored 25|, |26| as either an independent method to study air showers, or as a 
complement to existing fluorescence and surface detectors. A proposal has recently 
been put forth to evaluate the method using the Jicamarca radar system near Lima, 
Peru H- 

In order to increase the statistics at the high end of the spectrum significantly, two 
projects are now under preparation: 

(i) The Pierre Auger Observatory (PAO), currently under construction in Argentina, 
is the first experiment designed to work in a hybrid mode incorporating both a ground- 
based array of 1600 particle detectors spread over 3000 km^ with fluorescence telescopes 
placed on the boundaries of the surface array [2^. A second array will be set up in 
the Northern hemisphere to cover the whole sky. Such a full-sky coverage is very 
important to allow sensitive anisotropy analysis. The overall aperture (2 sites) for CRs 
with primary zenith angle < 60° and primary energy > 10^^ eV is ~ 1.4 x 10^ km^ sr. 

(ii) The mission "Extreme Universe Space Observatory" (EUSO) will observe the 
fluorescence signal of CRs, with energy > 4 x 10^^ eV, looking downward from the 
International Space Station to the dark side of the Earth atmosphere 29, li^. The 
characteristic wide angle optics of the instrument (with opening field of view ±30° at 
an average orbit altitude of ~ 400 km) yields a geometric aperture of ~ 5 x 10^ km^ 
sr, taking into account a 10% duty cycle. The monocular stand-alone configuration of 
the telescope will serve as a pathfinder mission to develop the required technology to 
observe the fluorescent trails of EASs from space. 

The experimental input for 7-rav physics will be further enriched by dedicated 



Cerenkov detectors like HESS 
as well as satellites like GLAST 



MAGIC tsa], CANGAROO |33|, and VERITAS [3 



35 



36 and AGILE [3 
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B. Primary species 



When a CR enters the Earth atmosphere it colhdes with a nucleus of an air atom, 
producing a roughly conical cascade of billions of elementary particles which reaches 
the ground in the form of a giant "saucer" traveling at nearly the speed of light. 

Unfortunately, because of the highly indirect method of measurement, extracting 
precise information from EASs has proved to be exceedingly difficult. The most fun- 
damental problem is that the first generations of particles in the cascade are subject to 
large inherent fiuctuations and consequently this limits the event-by-event energy reso- 
lution of the experiments. In addition, the center-of-mass energy of the first few cascade 
steps is well beyond any reached in collider experiments. Therefore, one needs to rely 
on hadronic interaction models that attempt to extrapolate, using different mixtures 
of theory and phenomenology, our understanding of particle physics. At present, the 
different approaches used to model the underlying physics of pp collisions show clear 
differences in multiplicity predictions which increase with rising energy 



Therefore, distinguishing between a proton and a nucleus shower is extremely difficult 
at the highest energies j41 



42 



Fortunately, photon and hadron primaries can be distinguished by comparing the 
rate of vertical to inclined showers, a technique which exploits the attenuation of the 
electromagnetic shower component for large slant depths. Comparing the predicted 
rate to the rate observed by Haverah Park for showers in the range 60° < 6 < 80°, 
Ave et al. ji^ conclude that above 10^^ eV, less than 48% of the primary CRs can be 
photons and above 4 x 10^^ eV less than 50% can be photons. Both of these statements 
are made at the 95% CL. 

The longitudinal development has a well defined maximum, usually referred to as 
-^max, which increases with primary energy as more cascade generations are required 
to degrade the secondary particle energies. Evaluating X^ax is a fundamental part of 
many of the composition studies done by detecting air showers. For showers of a given 
total energy, heavier nuclei have smaller X^ax because the shower is already subdivided 
into A nucleons when it enters the atmosphere. Specifically, the way the average depth 
of maximum (Xmax) changes with energy depends on the primary composition and 
particle interactions according to 

(Xmax) = -Deln ( ^ 



Eq 

where is the so-called "elongation rate" and Eq is a characteristic energy that 
depends on the primary composition ^]. Therefore, since (Xmax) and De can be 
determined directly from the longitudinal shower profiles measured with a fiuorescence 
detector, Eq and thus the composition, can be extracted after estimating E from the 
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total fluorescence yield. Indeed, the parameter often measured is Dio, the rate of 
change of (X^ax) per decade of energy. 

Another important observable which can be related to primary energy and chemical 
composition is the total number of muons N^^ reaching ground level. For vertical proton 
showers, numerical simulations j43] indicate that the muon production is related to the 
energy of the primary via 

E = 1.64 X 10^^ i^j eV. (2) 

Thus, modeling a shower produced by a nucleus with energy Ea as the collection of A 
proton showers, each with energy y4~^ of the nucleus energy, leads - using Eq. ((2)) - 

Consequently, one expects a CR nucleus to produce about 

^0.07 

more muons than a proton. This implies that an iron nucleus produces a shower 
with around 30% more muons than a proton shower of the same energy. 

The analysis of the elongation rate and the spread in X^ax at a given energy reported 
by the Fly's Eye Collaboration suggests a change from an iron dominated composition 



at 10^^'^ eV to a proton dominated composition near 10^^ eV j42l, l48| . Such behavior 
of De is in agreement with an earlier analysis from Haverah Park j49|. However, the 
variation of the density of muons with energy reported by the Akeno Collaboration 



51 



favours a composition that remains mixed over the 10^^ — 10^^ eV decade j^O]. More 

reanalyzed the Fly's Eye data considering not only 



recently, Wibig and Wolfendale 
proton and iron components (as in [4^) but a larger number of atomic mass hypotheses. 
Additionally, they adopted a different hadronic model that shifts the prediction of Xmax 
for primary protons of 10^^ eV from 730 g cm~^ [47] to 751 g cm~^. The difference, 
although apparently small, has a significant effect on the mass composition inferred 
from the data. The study indicates that at the highest energies (10^^'^ — 10^^ eV and 
somewhat above) there is a significant fraction of primaries with charge greater than 
unity. This result is more in accord with the conclusions of the Akeno group than 
those of the Fly's Eye group. Very recently, the Vo^ 
taking into account a bi-modal proton-iron model 
with 75 ± 5% of iron, with corresponding percentage of protons. A summary of the 
diferent bi-modal analyses is shown in Fig. ^ Within statistical errors and systematic 
uncertainties introduced by hadronic interaction models, the data seem to indicate 
that iron is the dominant component of CRs between ~ 10^^ eV and ~ 10^^ eV. 
Nonetheless, in view of the low statistics at the end of the spectrum and the wide 
variety of uncertainties in these experiments, one may conservatively say that this is 
not a closed issue. 



cano Ranch data was re-analyzed 
5^ 1. The best fit gives a mixture 
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FIG. 1: Predicted fraction of iron nuclei in the CR beam at the top of the atmosphere from 
various experiments: Fly's Eye (A), AGASA AlOO (■), AGASA Al (d) using sibyll 1.5 
as the hadronic interaction event generator ^] and Haverah Park [5J|, using QGSJET98 
(•) and QGSJEtOI (o) to process the hadronic collisions. The solid (dashed) line rectangle 
indicates the mean composition with the corresponding error estimated using the Volcano 
Ranch data and QGSJET98 (qgsjetOI); the systematic shift in the fraction of iron induced 
by the hadronic event generator is 14% 

C. Distribution of arrival directions 



The distribution of arrival directions is perhaps the most helpful observable in yield- 
ing clues about the CR orig in. On the one hand, if cosmic rays cluster within a small 
angular region (see e.g. [56[) or show directional alignment with powerful compact ob- 
jects (see e.g. j^), one might be able to associate them with isolated sources in the 
sky. On the other hand, if the distribution of arrival directions exhibits a large-scale 
anisotropy, this could indicate whether or not certain classes of sources are associated 
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with large-scale structures (such as the Galactic plane or the Galactic halo). 

Cosmic ray air shower detectors which experience stable operation over a period of a 
year or more can have a uniform exposure in right ascension, a. A traditional technique 
to search for large-scale anisotropics is then to fit the right ascension distribution of 
events to a sine wave with period 27r/m (m*^ harmonic) to determine the components 
(x, y) of the Rayleigh vector ji^ 

2 ^ 2 ^ 

x= — ^ cos(mai), V = sin(mai). (3) 

i=l i=l 

The m*^ harmonic amplitude of iV measurements of is given by the Rayleigh vector 
length Jl, = (x^ + y'^y^'^. The expected length of such a vector for values randomly 
sampled from a uniform phase distribution is TZq = 2/\fN. The chance probability 
of obtaining an amplitude with length larger than that measured is p{> TZ) = e~^°, 
where ko = 7^V7^o• To give a specific example, a vector of length ko > 6.6 would be 
required to claim an observation whose probability of arising from random fluctuation 
was 0.0013 (a "3(t" result) 3- 

AGASA has revealed a correlation of the arrival direction of the cosmic rays to 
the Galactic Plane (GP) at the Aa level js^. The energy bin width which gives the 
maximum fco-value corresponds to the region lO^'''^ eV - 10^^'^ eV where ko = 11.1, 
yielding a chance probability of p{> T^^^J^) ~ 1-5 x 10^^. The GP excess, which is 
roughly 4% of the diffuse flux, is mostly concentrated in the direction of the Cygnus 
region, with a second spot towards the Galactic Center (GC) [s^. Evidence at the 3.2cr 
level for GP enhancement in a similar energy range has also been reported by the HiRes 
Collaboration 6l|. The existence of a point-like excess in the direction of the GC has 
been confirmed via independent analysis of data collected with SUGAR.^ This is 
a remarkable level of agreement among experiments using a variety of techniques. 

At lower energies (~ PeV), the Rayleigh analysis shows no evidence of 
anisotropy j^]- Hence, the excess from the GP is very suggestive of neutrons as 
candidate primaries, because the directional signal requires relatively-stable neutral 
primaries, and time-dilated neutrons can reach the Earth from typical Galactic dis- 
tances when the neutron energy exceeds 10^^ eV. Arguably, if the Galactic messengers 
are neutrons, then those with energies below 10^^ eV will decay in flight, providing a 
flux of cosmic antineutrinos above 1 TeV that should be observable at kilometer-scale 
neutrino telescopes 65] . A measurement of the i/-fiux will supply a strong confirmation 
of the GP neutron hypothesis. 



^ Interestingly, sub- TeV 7-ray emission from the direction of the GC has been observed using the 
CANGAROO-II Imaging Atmospheric Cerenkov Telescope (6^. 
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For the ultra high en ergy 10 ^^'^ eV) regime, all experiments to date have reported 
ko -C 6.6, Vm < 5 0, |67|, |68|, |69|.^ This does not imply an isotropic distribution, but 
it merely means that available data are too sparse to claim a statistically significant 
measurement of anisotropy. In other words, there may exist anisotropics at a level too 
low to discern given existing statistics [tC^ . 

The right harmonic analyses are completely blind to intensity variations which de- 
pend only on declination, S. Combining anisotropy searches in a over a range of 
declinations could dilute the results, since significant but out of phase Rayleigh vectors 
from different declination bands can cancel each other out. Moreover, the analysis 
methods that consider distributions in one celestial coordinate, while integrating away 
the second, have proved to be potentially misleading [2l|. An unambiguous inter- 
pretation of anisotropy data requires two ingredients: exposure to the full celestial 
sphere and analysis in terms of both celestial coordinates. In this direction, a recent 
study [t^I of the angular power spectrum of the distribution of arrival directions of 
CRs with energy > 10^^'^ eV, as seen by the AGASA and SUGAR experiments, shows 
no departures from either homogeneity or isotropy on an angular scale greater than 
10°. Finally, the recently analyzed HiRes data is also statistically consistent with an 
isotropic distribution 



73 



All in all, the simplest interpretation of the existing data is that, beyond the ankle, 
a new population of extragalactic CRs emerges to dominate the more steeply falling 
Galactic population. Moreover, there are two extreme explanations for the near ob- 
served isotropy beyond 10^^'^ eV: one is to argue a cosmological origin for these events, 
and the other is that we have nearby sources (say, within the Local Supercluster) with 
a tangled magnetic field in the Galaxy, and beyond, which bends the particle orbits, 
camoufiaging the exact location of the sources. 

Although there seems to be a remarkable agreement among experiment on pre- 
dictions about isotropy on large scale structure, this is certainly not the case when 
considering the two-point correlation function on a small angular scale. The analyses 
carried out by AGASA Collaboration seem to indicate that the pairin g of events on the 
celestial sky could be occurring at higher than chance coincidence 56|,l2^- Specifically, 
when showers with separation angle less than the angular resolution 6*01111 = 2.5° are 
paired up, AGASA finds five doublets and one triplet among the 58 events reported 
with mean energy above 10^^'^ eV. The probability of observing these clusters by chance 
coincidence under an isotropic distribution was quoted as smaller than 1%. A third 



^ For the Fly's Eye data-sample the Rayleigh vector was computed using weighted showers, because 
it has had a nonuniform exposure in sideral time. A shower's weight depends on the hour of its 
sideral arrival time, and the 24 different weights are such that every time bin has the same weighted 
number of showers. 



11 



independent analysis using the Goldberg-Weiler formalism [t^, confirmed the re- 
sult reported by AGASA Collaboration and further showed that the chance probability 
is extremely sensitive to the angular binning. The "world" data set has also been stud- 
ied j?^. Six doublets and two triplets out of 92 events with energies > 10^^'^ eV were 
found, with the chance probability being less than 1% in the restricted region within 
±10° of the super-Galactic plane. The angular two-point correlation function of a com- 
bined data sarnEle of AGASA {E > 4.8 x 10^^ eV) and Yakutsk {E > 2.4 x 10^^ eV) 
was analyzed jlq]- For a uniform distribution of sources, the probability of chance 
clustering is reported to be as small as 4 x 10~^. Far from confirming what seemed a 
fascinating discovery, the recent analysis reported by the HiRes Collaboration showed 
that the data is consistent with no small-scale anisotropy among the highest energy 

IS- 



events 179 



The discovery of such clusters would be a tremendous breakthrough for the field, but 
the case for them is not yet proven. To calculate a meaningful statistical significance 
in such an analysis, it is important to define the search procedure a priori in order to 
ensure it is not inadvertently devised especially to suit the particular data set after 
having studied it. In the analyses carried out by AGASA Collaboration [H^E^, for 
instance, the angular bin size was not defined ahead of time. Very recently, with the aim 
to avoid accidental bias on the number of trials performed in selecting the angular bin, 
the original claim of AGASA Collaboration [s^ was re-examined considering only the 
events observed after the claim f8l]. This study showed that the evidence for clustering 
in the AGASA data set is weaker than was previously claimed, and consistent with the 
null hypothesis of isotropically distributed arrival directions. 

Summing up, the clustering on small angular scale at the upper end of the spec- 
trum remains an open question, and the increase in statistics and improved resolution 
attainable with PAO is awaited to solve the issue. 

D. Propagation of UHECRs 

In this section we briefly summarize the relevant interactions that CRs suffer on their 
trip to Earth. For a more detailed discussion the reader is refer to [l|, EE iSj, EB| • 

1. The GZK- cutoff 

Ever since the discovery of the cosmic microwave background (CMB) standard 
physics implies there would be a cuttofF in the observed CR-spectrum. In the mid- 



60's Greisen, Zatsepin, and Kuzmin (GZK) |86L l87l | pointed out that this photonic 



molasses makes the universe opaque to protons of sufficiently high energy, i.e., protons 
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with energies beyond the photopion production threshold, 

_ {nip + _ , ^ „ , , ^19 { £cmb 



E, 



th 

P7CMB 



6.8 X 10 



CMB 



10-3 eV 



-1 



eV 



(4) 



where nip (jriT,) denotes the proton (pion) mass and £^cmb ~ 10^'^ eV is a typical CMB 
photon energy. After pion production, the proton (or perhaps, instead, a neutron) 
emerges with at least 50% of the incoming energy. This implies that the nucleon en- 
ergy changes by an e-folding after a propagation distance < {ap^n^y)^^ ~ 15 Mpc. 
Here, n^ 410 cm-^ is the number density of the CMB photons, cTp^ > 0.1 mb is the 
photopion production cross section, and y is the average energy fraction (in the labo- 
ratory system) lost by a nucleon per interaction. Energy losses due to pair production 
become relevant below ~ 10^^ eV. For heavy nuclei, the giant dipole resonance can 
be excited at similar total energies and hence, for example, iron nuclei do not survive 
fragmentation over comparable distances. Additionally, the survival probability for 
extremely high energy {k. lO^'' eV) 7-rays (propagating on magnetic fields ^ 10^^^ G) 
to a distance ci, p{> d) ~ exp[— (i/6.6 Mpc], becomes less than 10"^ after traversing a 
distance of 50 Mpc.^ 

In recent years, several studies on the propagation of CRs (including both analytical 



analyses and nume rical simulations) have been carried out | 8^ iQOi . 



101 



102 
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104, 



10^ 



lEoi, 



107 



108 



95 



lOoT llol |. A summary of the 



|9£,^, 9a,M,iO 

UHECR attenuation lengths for the above mentioned processes (as derived in these 
analyses) is given in Fig. |21 It is easily seen that our horizon shrinks dramatically for 
energies > lO^'' eV. Therefore, if UHECRs originate at cosmological distances, the net 
effect of their interactions would yield a pile-up of particles around 4 — 5 x 10^^ eV with 
the spectrum droping sharply thereafter. As one can infer from Fig. the subtleties of 
the spectral shape depend on the nature of the primary species, yielding some ambiguity 
in the precise definition of the "GZK cutoff" . In what follows we consider an event to 
supersed the cutoff if the lower energy limit at the 95 % CL exceeds 7 x 10^^ eV. This 
conforms closely to the strong criteria outlined in Ref. jllll |. 



2. Propagation of CRs in a magnetized neighborhood of the Galaxy 

In addition to the interactions with the radiation fields permeating the universe, CRs 
suffer deflections on extragalactic and Galactic magnetic fields. 

Over the last few years, it has become evident that the observed near-isotropy of 
arrival directions can be easily explained if our Local Supercluster contains a large 



^ It should be stressed that if the extragalactic magnetic field is < 10 G, photons with enerj 
2> 10^^ eV can reach us without significant energy loss from distant (redshift z > 0.03) sources [8] 
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FIG. 2: Attenuation length of 7's, p's, and ^^Fe's in various background radiations as a 
function of energy. The 3 lowest and left-most thin solid curves refer to 7-rays, showing the 
attenuation by infra-red, microwave, and radio backgrounds. The upper, right-most thick 
solid curves refer to propagation of protons in the CMB, showing separately the effect of 
pair production and photopion production. The dashed-dotted line indicates the adiabatic 
fractional energy loss at the present cosmological epoch (see e.g. Appendix B of Ref. 
The dashed curve illustrates the attenuation of iron nuclei. 



scale magnetic field which provides sufficient bending to the CR trajectories lid . 113 
Intergalactic field strengths and coherence lengths are not well established, but it is 
plausible to assume that fields have coherent directions on scales i ~ 0.5 — 1 Mpc. The 
Larmor radius of a CR of charge Ze propagating in a magnetic field B^g = G 
is given by 

100 E20 



Z B 



Mpc 



nG 



(5) 
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where £"20 is the particle's energy in units of 10^° eV. For tl ^ £ the motion is not very 
different from a quasihnear trajectory, with small deflections away from the straight 
line path given by 

^ Q 30 W^A^ ^ (6) 



E20 

where Lmpc is the propagation distance in units of Mpc. As the Larmor radius starts 
approaching i the particles begin to diffuse. 

Diffusion has two distinctive regimes. Particles that are trapped inside magnetic 
subdomains (of size ^Mpc = ^/Mpc) follow Kolmogorov diffusion. In such a ca se, the 
functional dependence of energy of the difussion coefficient is found to be |ll4l | 



D{E) ^ 0.0048 



-^20 ^Mpc 



1/3 



MpcVMyr . 



(7) 



With rising energy, tl 



and there is a transition to Bohm diffusion. The diffusion 



1151. In 



coefficient in this regime is of order the Larmor radius times velocity (~ c) 
this case the accumulated deflection angle from the direction of the source, ca n be 
estimated a.=uming that the particle, make a random walk in the magnetic field Ei 



e(E) » 0.54= 



f/1/2 J 1/2 y ry 

'Mpc -^Mpc ^ -"nG 



E. 



(8) 
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Surprisingly little is actually known about the extragalactic magnetic field strength. 
There are some measurements of diffuse radio emission from the bridge area between 
Coma and Abell superclusters that under assumptions of equi parti tion allows an es- 
timate of 0.2 — 0.6 /iG for the magnetic field in this region [ll7|.^ Such a strong 
magnetic field (wh ich is compatible with existing upper limits on Faraday rotation 
measurements |ll9|) could be possibly understood if the bridge region lies alonga fila- 



ment or sheet of large scale structures |l2ol |. Faraday rotation measurements |lig 
have thus far served to set u pper boun ds of O{10^^ — 10^^) G on extragalactic mag- 
netic fields on various scales lll[ 
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as have the limits on distortion of the CMB 
The Faraday rotation measurements sample extragalactic field strengths of 
any origin out to quasar distances, while the CMB analyses set limits on primordial 
magnetic fields. Finally, there are some hints suggesting that the extrag alact ic field 



strength can be increased in the neighborhood of the Milky Way, B^q > 10 125l |. Now, 
using Eq. one can easily see that because of the large uncertainty on the magnetic 
field strength, C(nG) — C(/iG), all 3 different regimes discussed above are likely to 
describe UHECR propagation. 



Fields of 0{fj.G) are also indicated in a more extensive study of 16 low redshift clusters 
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If CRs propagate diffusively, the radius of tlie spliere for potential proton sources 
becomes significantly reduced. This is because one expects negligible contribution 
to the flux from times prior to the arrival time of the diffusion front, and so the 
average time delay in the low energy region, rjeiay ~ d'^ /[^D{E)], must be smaller than 
the age of the source, or else the age of the universe (if no source within the GZK 
radius is active today, but such sources have been active in the past). Note that the 
diffuse propagation of UHE protons requires magnetic fields ~ l/^G. Therefore, for 
typical coherence lengths of extragalactic magnetic fields the time delay of CRs with 
E ^ 10^^'^ eV cannot exceed r^eiay ^ 14 Gyr, yielding a radius of (i ~ 30 Mpc. In 
the case CR sources are active today, the radius for potential sources is even smaller 

~ 5 Mpc. 

On the other hand, the sphere of potential nucleus-emitting-sources is severely con- 
strained by the GZK cutoff: straightforward calculation, using the attenuation length 
given in Fig.|21 shows that less than 1% of iron nuclei (or any surviving fragment of their 
spallations) can survive more than 3 x 10^^ s with an energy > lO^'^'^ eV. Therefore, 
the assumption that UHECRs are heavy nuclei implies ordered extragalactic magnetic 
fields 5nG ^ 15 — 20, or else nuclei would be trapped inside magnetic subdomains 
suffering catastrophic spallations. 

The large scale structure of the Galactic magnetic field carries substantial un- 
certainties as well, because the position of the solar system does not allow global 
measurements. The average field strength can be directly determined from pul- 
sar observations of the rotation and dispersion measures average along the line of 
sight to the pulsar with a we igh proportional to the local free electron density, 
f» 2/iG 12fil . 127 . 128 . 129^ . (We use the standard, though ambiguous notation, in 



which B refers to either the Galactic or extragalactic magnetic field, depending on the 
context.) Measurements of polarized synchrotron radiation as well as Faraday rotation 
of the radiation emitted from pulsars and extragalactic radio sources revealed that the 
global structure of the magnetic field in the disk of our Galaxy could be well desc ribed 
by spiral fields with 27r (axisymmetric, ASS) or vr (bisymmetric, BSS) symmetry 130l . 



In the direction perpendicular to the Galactic plane the fields are either symmetric 
(S) or antisymmetric (A). Discrimination between these models is complicated. Field 
reversals are certainly observed (in the Crux-Scutum arm at 5.5 kpc from the Galactic 
center, the Carir ia-Sa gittarius arm at 6.5 kpc, the Perse us a rm at 10 kpc, and possibly 
another beyond |l3lj ). However, as discussed by Vallee |l32| . turbulent dynamo theory 
can explain field reversals at distances up to ~ 15 kpc within the ASS configuration. 
More accurately, the fi eld streng th in the Galactic plane {z = 0) for the ASS model 



is generally described by |l33l Il34 



Bip,e) = Boip) cos2[0-/51n(p/eo)], (9) 
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and for the BSS 

B{p,e) = B,{p) cos[^-/31n(p/eo)], (10) 

where 9 is the azimuthal coordinate around the Galactic center (clockwise as seen from 
the north Galactic pole), p is the galactocentric radial cylindrical coordinate, and 

5o(p) = — tanh3(p/pi) pG. (11) 
P 

Here, = 10.55 kpc stands for the galactocentric distance of the maximum of the field 
in our spiral arm, (3=1/ tanp (with the pitch angle, p = —10°), tq = 8.5 kpc is the 
Sun's distance to the Galactic center, and pi = 2 kpc. The 6 and p coordinates of the 
field are correspondingly. 

Be = B{p, 9) cos p , B^ = B{p,9) sinp . (12) 

The field strength above and below the Galactic plane (i.e., the dependence on z) has 
a contribution coming from the disk and another from the halo: (i) for A models 



5a (p, 9, z) = Bip, 9) tanh(z/2;3) — — - + r-— - , (13) 

\2cosh(z/zi) 2 cos\i[z I Z2) J 

(a) for S models, i?s = Bx{p,9, z) /taiih.{z/ z^); where zi = 0.3 kpc, ^2 = 4 kpc and 
z^ = 20 pc. With this in mind, th e Ga lactic magnetic field produce significant bending 
to the CR orbits if E20/Z = 0.03 
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E. GZK-end of the cosmic ray spectrum? 



A first hint of a puzzle surfaced in the highest energy F ly's Eye event |135| which 
has no apparent progenitor within the L ocal Supercluster [l36[. Subsequent observa- 
tions with the AGASA experiment |l37| carried strong indication that the cutoff was 
somehow circumvented in the absence of plausible nearby sources. 

T he big disa ppointment of 2002 was the CR-flux reported by th e HiR es Collabora- 



tion 



138L 



139| . which is in sharp disagreement with AGASA data 14ol |. The discrep- 
ancy between the two estimated fluxes is shown in Fig. El One can argue correctly 
that the statistical significance of the discrepancy is small, although such an assesment 
requires a conspiracy between the two groups to be nd their maximal systematic errors 
in opposite directions. Moreover, an analysis 143| of the combined data reported by 
the HiRes, the Fly's Eye, and the Yakutsk collaborations is supportive of the existence 
of the GZK cutoff at the > 5a (> 3.7a) level. The deviation from GZK depends on 
the set of data used as a basis for power law extrap olati on from lower energies. An 
additional input for this analysis was the recent claim |l44l | that there may be technical 



17 



26 



25,5 



25 



O 



■ SUGAR 

• AGASA 

A Fly's Eye (stereo) 

* Haverah Park 
VHiRes— 1 (monocular) 



24.5 



24 



1 i 



{ 



1 



23,5 



23 



19 19,2 19,4 19.6 19,8 20 20,2 20.4 20.6 

Log E (eV ) 

FI G. 3: Upper end of t he co smic ray energy spectrum as observed by AGASA fl4^, Fly'; 
Eye|l4l|, Haverah Park [uil, HiRes and SUGAR 0- 



problems with the Yakutsk data collection. More recently, fingerprints of super-GZK 
CRs have been found [i^ by reanalyzing the SUGAR data However, as one can 
see in Fig. |3J the number of events is not enough to weight in on one side or the other 
with respect to the GZK question. 



II. VANILLA SKY: UHECR GENERATION WITHIN THE STANDARD 
LORE 

A. A brief low energy perspective 

Supernova remnants (SNRs) are thought to be the main source of both CR ions 
and electrons with energies below the knee. The particle acceleration mechanism in 
individual SNRs is usually assumed to be diffusive shock acceleration, which naturally 
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leads to a power-law popu lation of relativistic particles. In the standard version of 



this mechanism (e.g. 145|), particles are scattered by magnetohydrodynamic waves 



repeatedly through the shock front. If they encounter an enhancement of mole cular 
density, the pion channel can lead to observable amounts of 7-rays (see Ref. 
for a review, and references therein for details). Electrons suffer synchrotron losses, 
producing the non-thermal emission from radio to X-rays usually seen in shell-type 
SNRs. The maximum energy achieved depends on the shock speed and age as well as 
on any competing loss processes. In young SNRs, electrons can easily reach energies 
in excess of 1 TeV , where they produce X-rays by synchrotron mechanism (see, for 
example, |i47lE48| 1. 

CRs of low energies are also expected to be accelerated in OB ass ociations, through 
turbulent motions and collective effects of star winds (e.g. 149 . 150l |). The main accel- 
eration region for TeV particles would be in the outer boundary of the supperbubble 
produced by the core of a given stellar association. If there is a subgroup of stars 
located at the acceleration region, their winds might be illuminated by the locally ac- 
celerated protons, which would have a distri butio n with a slope close to the canonical 
value, a ~ 2, and produce detectable T- ravs 15 1| . The HEGRA detection in the vicin- 
ity of Cygnus 0B2, TeV J2032-t-4131 152j |. could be, judging from multiwavelength 
otee™aons|I3, the result of =uch a proc ess Il5l| . Anearby EGRET source (3EG 
J2033+4118) has also a likely stellar origin |l54 lisd [l5fi| . 

Truth is, as always, bitter. No astrophysical source of UHECRs, nor of CRs with 
energies below the knee, has been ever confi rmed. Out of all SN Rs coinciding with non 



variable 7-ray sources detected by EGRET |l53, EHl EHl H , 

the supernova remnant 



RX J1713. 7-3946 is perhaps one the most convincing cases for a hadronic cosmic- 
rav a ccele r ator dete c ted so far in the Galaxy, although yet subject to confir matio n 



I6IL II62L I163L I164L 16^ . Other excellent candidates include SN1006 (e.g. 



and Gas A (e.g. |l67 |). Kilometer-scale neutrino t elescopes have also been proposed 



168, 



169, 07 
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and will be 



as viable detectors of hadronic GR sources (e.g. 
a welcomed addition to the space- and ground-based detectors already existing or 
planned. 



B. Plausible sources of UHECRs and the Hillas' plot 

Following the main ideas behind the concept of Fermi's first order acceleration, when 
tl approaches the accelerator size, it becomes very difficult to magnetically confine the 
GR to the acceleration region, and thus to continue the accelerating process up to higher 
energies. If one includes the effect of the characteristic velocity f3c of the magnetic 
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FIG. 4: The Hillas diagram showing (chain curves) magnetic field versus gyroradius for 
proton momenta 10^^, 10^^, . . . , 10^^ eV/c. The solid curves correspond to different shock- 
waves velocities: the upper solid curve indicates the maximum attainable energy (3 = 1, the 
middle and lower solid curves indicate plausible less effective acceleration processes. Typical 
size and magnetic field of possible acceleration sites are shown for neutron stars (ns), white 
dwarfs (wd), sunspots (ss), magnetic stars (ms), active galactic nuclei (ag), interstellar space 
(is), SNRs (sn), radio galaxy lobes (rg), galactic disk (d) and halo (h), clusters of galaxies 
(cl) and intergalactic medium (ig). Typical jet-f rame parameters of the synchrotron proton 
blazar model 



174l | and gamma ray burst model '1251 are indicated by open squares labeled 



"bl" and "gb" , respectively [1 
scattering centers^, the abo ve a rgument leads to the general condition (sometimes 



called the "Hillas criterion" [122]) 



E^ax~2/3cZe5r,, (14) 

for the maximum energy acquired by a particle travelling in a medium with magnetic 
field B. 



^ The size of the accelerating region containing the magnetic field should be as large as 2r^^ . Taking 
into account a characteristic velocity (3c of the scattering centers this transforms into 2t^I(3. 
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In the case of one-shot acceleration scenarios, the maximum reachable energy turns 
out to have a quite similar expression to the shock acceleration case of Eq. (fT^. For in- 
st dimensional analysis suggests that the maximum energy that can be obtained 

from a pulsar is 
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E^^^ = -ZeB,Rl, (15) 
c 

where uj is the pulsar angular velocity, Bg the surface magnetic field and Rg the neutron 
star radius. Therefore, if Eg ~ 10^^ G, Rg ~ 10 km, and uj ~ GOvr s~^ (as for the Crab 
pulsar), a circuit connected between pole and equator would see an emf ~ 10^*^ V for an 
aligned or oblique dipole. When realistic models of acceleration are constructed, how- 
ever, this ideal dimensional limit is not fully realized, because the large potential drop 
along the magnetic field lines is significantly short-circuited by electron and positrons 
moving in opposite directions along the field lines |l73l |. 

The dimensiona l arg uments of Eqs. (fT^ and ()15|1 are conveniently depicted in the 



"Hillas diagram" 172l | shown in Fig. 0] Very few sites can generate particles with 
energies > lO^'' eV: either this occurs on highly condensed objects with huge B or 
enormously extended objects. Some of these potential astrophysical sources are dis- 
cussed in what f ollows. F or further details on the electrodynamical limitations of CR 



sources see, e.g. 



177, 17 



C. Neutron stars 

1. Magnetohydronamic acceleration of iron nuclei in pulsars 

Following earlier ideas |l79l |^. Blasi et al. |l82l | ^ have shown that young magnetized 
neutron stars in our own Galaxy may be one such astrophysical system that satisfies 
the Hillas' criterion. 

Neutron stars -endpoints of stellar evolution- begin their life rotating rapidly {Q ~ 
3000 rad s~^) and with large surface magnetic fields (Bs > 10^^ G). The dipole compo- 
nent decreases as the cube of the distance from the star's surface B{r) = Bs{Rs/rY, 
where the radius of the star is Rs — L. At the light cylinder Ric = c/^7 ~ l^'^Vt'^^ cm, 
where Vt^^k = i^/SOOO rad s~^, the dipole field cannot b e cas ually maintained, the field 
is mostly azimuthal, with field lines spiraling outwards |l87| . Inside the light cylinder. 



^ Even earlier ideas relating UHECRs with neutron stars can be found in Refs. |l8d 1X8 1| . although 
these attempts have failed at either reaching the highest energies, or reproducing the spectrum, or 
repr pduf-iTig t he apnaient isotropy of the arrival di rect i nns of UHECRs. 

^ See |l83lll84l| . also |l85| . for related research. See |l86l | for yet another model of UHECR generation 
in neutron stars, involving planetoid impacts. 
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the magnetosphere corotates with the star, and the density of material (mostly iron 
peak elements formed during the supernova event that were stripped off the surface 
due to strong electric fields) has t he Goldreich- Julian value, ncjir) = B{r)Q/ {inZec), 
where c is the speed of light llSS^ . The j p ehayior of the plasma outside the light cylin- 
der is still not yet fully understood isd IQfll 191 . 192| . although some analytical and 



numerical studies show the development of kinetically dominated relativistic winds (see 
e.g., H). 



Blasi et al. 1821 ] assumed that the magnetic field in the wind zone decreases as 
B{r) < BicRic/r. For surface fields of Bs = 10^^ i?i3 G, the field at the light cylinder 
is Bic = 10^° Bi^Q^i^ G. The maximum energy of particles that can be contained in the 
wind near the light cylinder is 

i^max = ^ 8 X 10^° Z2eB,,nl, eV , (16) 

c 

where = Z/26. The typical energy of the accelerated CRs, E^r, can be estimated 
by considering the magnetic energy per ion at the light cylinder Ecr — Bf^/87mQj. At 
the light cylinder ncj = 1-7 x 10^^ Bi^Q'^f^/Z cm""^ which gives 

E„ ~ 4 X 10^0 Z2eBi3Qlk eV . (17) 

In this model, as the star spins down, the energy of the cosmic ray particles ejected 
with the wind decreases. The total fiuence of UHECRs between energy E and E + dE 
is 

N{E)dE = ^^dE, (18) 
where the particle luminosity is 

Ar = CnGj TiRlc = 6 X lO^^^^lfk (19) 

Z26 

and ^ < 1 is the efficiency for accelerating particles at the light cylinder. For a spin 
down rate dominated by magnetic dipole radiation, given by IQQ = —BgRgQ^/6c^ 
where I = 10^^ g cm^ is the moment of inertia, the time derivative of the spin frequency 
ish = 1.7 X 10-^523^3^ s-S and Eq. (HH) gives 

Substituting in Eq. ()18p. the particle spectrum from each neutron star is 

N{E) = d^^T GeV' . (21) 
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Taking the confining volume for tliese particles to be Vc and the lifetime for con- 
finement to be tc, the UHECR density is n{E) = eN{E)tc/TVc, and the flux at the 
surface of the Earth is F{E) = n{E)c/4:. For a characteristic confinement dimension 
of i? = 10 Ri kpc, Vc = iivR^/S and = QR/c, where Q > 1 is a measure of the how 
well the UHE CR are trapped. With this in mind, the predicted UHECR flux at the 
Earth is [isi ^ 



F{E) = 10-'^ GeV-^cm-V^ • (22) 

Here, the fact that neutron stars are produced in our Galaxy at a rate 1/r, where r = 
100 T2 yr, but that not all them (but rather only a fraction e) have the required magnetic 
fields, initial spin rates and magnetic field geometry to allow efficient conversion of 
magnetic energy into kinetic energy of the flow, was taken into account. By comparing 
with observations, the required efficiency factor, ^e, can be estimated, and it only needs 
tobe^e>4 X 10-^Q-\ 

The condition that a young neutron star could produce the UHECRs^ is that Ecr 
exceeds the needed energy when the envelope becomes transparent (i.e. before the 
spinning rate of the neutron star decreases to the level where the star is unable to 
emit particles of the nec essary energy), EcriUr) > 10^°i?2o eV. This translates into the 
following condition |l82l | 

^. > — ^ . (23) 



r1/2 
-Dl3 



4^26^20' -0.13Mi5i3^5i'/' 



Eq. (j23j) translates in turn into upper bounds on the surface magnetic field strength 
and the star initial spin period Pi = 27r/cijj, 

i?i3 < ^i^H^ , (24) 



® Note that the predicted spectrum of Eq. (|22|l is very flat, 7=1, which would agree only with 
the lower end of the plausible range of 7 observed at ultra-high energies. Propagation effects can 
produce an energy dependence of the confinement parameter Q and, correspondingly, a steepening 
of the spectrum toward the middle of the observed range 1 ^ 7 ^ 2. 

^ The gyroradius of these UHECRs in a Galactic field of /iG strength is considerably less than the 
typical distance to a young neutron star 8 kpc). Since the iron arrival distribution at 10^° eV 
probes similar trajectories to protons at a few times 10^^ eV, Galactic iron nuclei would show a 
nearly isotropic distribution with a slight correlation with the Galactic center and disk, at higher 
energies. 

° A supernova that imparts Esn = 10^-^ ^51 erg to the stellar envelope of mass Menv = 10 Mi M0 
is considered. Also, that the condition fo r iro n nuclei to traverse the supernova envelope without 
significant losses is that E < 100 g cm~^ |l93| . 
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and 

P,<8nBli^Z^,E^oK (25) 

For Ml = 2 and E20 = E^i = Z26 = 1, Eq (glj) gives B13 < 15.4, whereas Eq. ^ leads 
to Pi < 10 ms, not very restrictive values for a young neutron star, see for example the 
Parkes Multibeam Pulsar Survey. 



2. Magnetars 

Magnet ars , neutron stars with surface dipole fields on the order of 10^^ G llQ'i 

l20l| , were also proposed as plausible sites for the generation of UHECRs [20^ 



Assuming that they occur in all galaxies which form massive stars (then avoiding the 
large-distance GZK problem), and that the UHECR are arriving from outside our own 
Galaxy, the lumin ous infrared galaxies are preferred sites to search for a magnetar 
origin of CRs (see 



2031 1 and below the discussion in Section Uli 

IS 
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The magnetar model for the acceleration of UHECRs proposed by Arons 
a variant of that using neutron stars outlined in the previous Section. The theory 
predicts an injection of charged particles with maximum energy 

E^ax = Ze^i = Ze^^ = 3 x lO'^^ZB^^nl eV, (26) 

where i?* is a magnetar's surface magnetic field, B15 = i?^,/(10^^ G), Qi = 10'^f24 s"^ 
is the initial angular velocity of the neutron star, is the stellar radius, and c is the 
speed of light. The initial rotation period is Pi = O.6A/Q4 ms {if Z = 1 — 2, one requires 
Pi < 2 — 3 ms for the model to be viable). 

The ions actually gain their energy in the relativistic wind electromagnetically ex- 
pelled from the neutron star at distances r larger than the radii of the star and its 
magnetosphere. This avoids catastrophic radiation losses; the electric potential in the 
wind is rE = rB = $. As the star spins down, as in the model by Blasi et al., 
the voltage and the maximum particle energies decline. Summing over the forma- 
tion and spindown event, one finds a per event injection spectrum proportional to 
f{E) = E^^[l + (E/Eg)]^^ for E < -Emax- Here Eg measures the importance of grav- 
itational wave losses (calculated for a star with static non-axisymmetric quadrupole 



The Parkes multibeam pulsar survey is a large-scale survey of a narrow strip of the inner Galactic 
plane (|6| < 5°, 260° < I < 50°, see and references therein). It has much greater sensitivity 



than any previous survey to young and distant pulsars along the Galactic plane, and it has resulted 
in the detection of many prey ii 
7-ray sources, e.g. |l95lll96| . 



in the detection of many previously unknown young pulsars, potentially counterparts of unidentified 
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asymmetry) in spinning the star down. When they exert torques larger than the elec- 
tromagnetic torque, the star spends less time at the fastest rotation rates (i.e. less 
time accelerating the highest energy particles), thus causing a steepening in the spec- 
tral slope at the highest energies. If the star has an internal magnetic field even stronger 
than the already large surface field, equatorial ellipticities in excess of 10~^ can exist, 
in which case would be less than -Emax- Three cases of energy loss due to gravi- 
tational radiation (GR) were considered in the model: no GR loss (ee = 0, £'g = oo); 
moderate GR loss (ee = 0.01, Eg = 3 x lO^o eV); strong GR loss (ee = 0.1, Eg = 3 x 10^^ 
eV). If one assumes that all magnetars have exactly the same starting voltage (10^^'^ V), 
then the model predicts that the spectrum E^J{E) should rise with E above the en- 
ergy Eg = 2.8 X 10^° eV, where the GZK loss rate becornes approximately energy 
independent (unless the gravitational wave losses are large) |202| . This prediction will 
certainly be testable with the Pierre Auger Observatory. 

The total number of particles injected per event is 

c^R^I 10^3 
^ ^ (27) 

for a stellar radius of 10 km and a moment of inertia / = 10*^^ cgs. The rate at which 
galaxies inject UHECR into t he u niverse in this model then is ricr = ^'mag^i'^gaiaxy, 
where ngaiaxy ~ 0.02 Mpc"^ |2o4 . iVi is given by Eq. (j2Z|), and z/^fg is the birth 
rate of rapidly rotating magnetars per galaxy. Multiplication of the source spectrum 
q{E) oc hcrf{E) by the energy dependent GZK loss time yields a spectrum received 
at the Ear th in reasonable accord with the existing observations of UHECR if z/J?^* ^ 

I I mag 

10^^ yr~^ (2021 . That fast magnetar birth rate lies between 1% and 10% of the total 
magnetar birth rate inferred for our galaxy, and about 0.1% of the total core collapse 
supernova rate in average star forming galaxy, ~ 10~^ yr~^ flO^. 



3. UHECRs from a pulsar in Cygnus 0B2? 



As discussed in Sec. Ill Al some evidence may be emerging for a GR accelerator in 
the Gygnus spiral arm. The HEGRA experiment has detected an extended TeV 7- 
ray source in the Gygnus region with no clear counterpart and a spectrum not easily 



accommodated with leptonic radiation 152| . The difficulty in accommodating the 



spectrum by conventional electromagnetic mechanisms has been exacerbated by the 
failure of Ghandra and VLA t o det ect significant levels of X-rays or radiowaves signaling 
acceleration of any electrons 1153^ . Especia lly i ntriguing is the possible association of 
this source with part of Gygnus 0B2 itself |l5l|. a cluster of several thousands young, 
hot OB stars with a total mass of ~ 10^ j206[. At a relatively small distance to 
Earth, f» 1.7 kpc, this is the largest massive Galactic association. It has a diameter of 
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~ 60 pc and a core radius of ~ 10 pc. The typical main sequence evolution lifetime of 
massive O stars is O (Myr) and a few tens Myr for massive B stars. Since the 0-star 
population should pass through the Wolf-Rayet phase and explode as supernovae, very 
fast pulsars are expected to be born in explosions of these massive stars at a rate of 
about one every ten thousand years. 

Apart from the mentioned interpretation of a hadronic production of the TeV ra- 



diation within the winds of outlying OB stars of Cyg 0B2 |151l |. it was recently put 
forward that the TeV emission reported by HEGRA and the CR anisotropy observed 
at about 10^^ eV in the direction of the Cygnus region can be related to a young pulsar 
and its pulsar wind nebu lae (PWN), born in the Cygnus 0B2 association a few ten 
thousands years ago |207l |. The TeV 7-ray emission would originate in the PWN as a 
result of interactions of high energy hadrons and/or leptons, whereas there would be 
a directional CR signal due to neutrons that are dissolved from heavy nuclei acceler- 
ated by the pulsar. Within this model, however, it is hard to explain the absence of 
counterparts at lower (EGRET) energies at the location of the TeV source, as well as 
the absence of a stronger X-ray source. However, disregarding these difficulties, an 
interesting consequence for CR physics can be inferred. 

As already discussed in Sec. Ill C 1\ heavy nuclei can attain ultra high energies through 
magnetic sling shots inside the neutron star wind zone. The small part of nuclei, which 
escaped from the pulsar wind nebula, are likely to be captured by strong magnetic fields 
of dense regions of the OB associat ion. Magnetic field strengths in dense molecular 



clouds are expected to b e ~ 1 mG [209|. Thus, nuclei with E/Z 1 EeV prop agate 



attaining Bohm diffusion |207l | . The resulting time delay of several thousand years 



produces a steepening of the characteristic power law injection spectrum, oc 
of the pulsar nebula. In their random traversal of the OB association, the nuclei 
undergo photodisintegration on the far infrared thermal photon population and liberate 
neutrons. 



A similar scenario would explain the anisotropy spot towards the Galactic Center 2Q8J. 
Indeed, location problems may arise as well: the PWN size (given that the TeV source size is 6 
pc if indeed located at 1.7 kpc, and that the source is diffuse) would make the PWN substantially 
larger than both Vela's PWN (~0.1 pc at 250 pc) and Crab's (~1 pc). Additionally, if one were to 
associate the nearby EGRET source 3EG 2033+4118 with the putative pulsar itself, it is unclear 
whether the PWN hypothesis for the TeV source would imply that it is only a one sided PWN (Y. 
Butt, private comunication) . 

Note that once diffusion has been established, additional Rayleigh steps in the Galactic magnetic 
field do not change the spectral index —3) significantly. 
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The thermal photon density at the source reads, 



T3 x^dx (4.37 r)3 x^dx , 
'^IR = — / T ~ ^ / 7 cm -\ {21 



where T is the kinetic temperature of the molecular cloud, x = E^/T, and E^ is the 
target photon energy in degrees K. The relevant photon energy range is established via 
standard kinematics of the photodisintegration process. 



i5, = §^ =10|^K. (29) 



where 5 < -EmcV < 25 is the energy region for giant dipole resonance contribution 
in the nucleus rest frame, and E^/A = Ejsf Pey is the lab frame energy /nucleon. For 
EN,PeY ~ 10^ Eq. (121 leads to Xi = 50/T and xs = 250/T, with T in degree Kelvin. 
Molecular clouds with HII regions have temperatures between 15 and 100 K 21l|, 
thus taking an average photodisintegration cross section of 40 mb,^^ we find that the 
nucleus mean free path lies between 0.80 and 380 pc. This corresponds to a reaction 
time between 4 to 1500 yr -C time delay, allowing sufficient neutron production to 
explain the anisotropy |65l |. 

The Galactic anisotropy observed by the various collaborations spans the energy 
range 0.8 to 2.0 EeV. The lower cutoff specifies that only neutrons with EeV energies 
and above have a boosted cr„ sufficiently large to serve as Galactic messengers. The 
decay mean free path of a neutron is cr„r„ = 10 (i?„/EeV) kpc, the lifetime being 
boosted from its rest-frame value r.„ = 886 seconds to its lab value via r„ = i?„/m„. 
Actually, the broad scale anisotropy from the direction of the GP reported by Fly's Eye 
Collaboration peaks in the energy bin 0.4 — 1.0 EeV, but persists with statistical 
significance to energies as low as 0.2 EeV. This implies that if neutrons are the carriers 
of the anisotropy, there needs to be some contribution from at least one source closer 
than 3-4 kpc. Interestingly, the full Fly's Eye data include a directiona l signal fr om 
the Cygnus region which was somewhat lost in unsuccessful attempts |212l . l213| to 
relate it to 7-ray emission from Cygnus X-3. The upper cutoff refiects an important 
feature of photodisintegration at the source: heavy nuclei with energies in the vicinity 
of the ankle will fragment to neutrons with energies about an order of magnitude 
smaller. To account for the largest neutron energies, it may be necessary to populate 
the heavier nucleus spectrum in the region above the ankle. This is not a problem 
- one fully expects the emerging harder extragalactic spectrum to overtake and hide 



The photoabsorpti on cr oss section roughly obeys the Thomas- Reiche-Kuhn sume rule, i.e., = 
60NZ/A mb-MeV [u^. 

To produce the highest energy neutrons (with energy < 10^*'^ eV) via photodisintegration of 
medium mass (say, A = 10 — 20) nuclei, one needs primary particles with energies < 10^^'^ eV. For 
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the steeply falling galactic population. It is not therefore surprising that in order to 
fit the spectrum in the anisotropy region and maintain continuity to the ankle region 
without introducing a cutoff, the AGASA Collaboration required a spectrum oc 
or steeper ji^. 

For every surviving neutron at ~ EeV, there are many neutrons at lower energy that 
decay via n ^ p + e~ + V^. The proton is bent by the Galactic magnetic field, the 
electron quickly loses energy via synchrotron radiation, and the ^travels along the 
initial neutron direction, producing a directed TeV energy beam j^. The sensitivity 
of forthcoming neutrino telescopes to this signal is discussed in Sec. IIVBI 



D. Radio Galaxies and Active Galactic Nuclei 

1. Definitions 

Blazars are active galactic nuclei (AGNs) with a) strong flat spectrum radio emission 
[the power law index a > —0.5, with S'(z/) oc z/"] and/or b) significant optical polar- 
ization, and/or c) significant flux variability in the optical and in other wavelengths. 
When the optical variability occurs on short timescales, the objects are referred to 
as optically violently variable -OVV- quasars. The blazar classification also includes 
BL Lacertae (BL Lac) objects, which present a complete or nearly complete lack of 
emission lines, and highly polarized quasars (HPQs). It also refers, sometimes, to fiat 
spectrum radio quasars (FSRQs), although these are generally more distant, more lu- 
minous, and have stronger emission lines. Within the unification model, the underlying 
scenario for all AGNs is intrinsically similar. At the very center of the galaxy there is 
a supermassive black hole (~10® to ~10^° M©) which accretes galactic matter forming 
an accretion disk. Broad emission lines are produced in clouds orbiting above the disc 
at high velocity, the broad line region (BLR) and this central region is surrounded by 
an extended, dusty, molecular torus. A hot electron corona populates the inner region, 
probably generating continuum X-ray emission. Narrower emission lines are produced 
in clouds moving much farther from the central black hole. Two-sided jets of relativis- 
tic particles emanate perpendicular to the plane of the accretion disc, the generation 
of which is still not fully understood. Unification of different AGN classes is achieved 



a falling spectrum cx E^^, the medium mass nucleus population at 10^^'^ eV is roughly 3 orders 
of magnitude smaller than the diffuse flux at 10^^'^ eV. This means that it is about 1.6 orders of 
magnitude smaller than the reported CR excess (which is about 4% of the diffuse flux, see Sec. II C|l . 
However, since each nucleus produces roughly A — Z neutrons, the average number of liberated neu- 
trons is on the same order of magnitude, 1.6 — log]^Q(A — Z), than the obseved neutron population 
at < 101^-2 eV. 
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FIG. 5: The unification model for AGNs. The components of the figure are discussed in the 
text. Blazars are those AGNs for which the jets are close to line of sight. A regular quasar 
or a Seyfert 1 galaxy is observed if the orientation angle is ~ 30°, where the narrow- line 
and broad-line regions are visible. At larger angular offsets, the broad-line region will be 
hidden by the torus, the corresponding class being Seyfert 2 galaxies. Perpendicular to the 
jet axis, the full extent of the jets may be seen particular at low freque ncies , giving rise to a 
morphology typical of radio galaxies. The figure is adapted from Refs. Q, Ei, Ei, Q . 



taking into account t he in trinsic anisotropy of the phenomenon, is shown in Fig. El (see 



214 . 215 . 216L 217 1 for further and more detailed discussions) 



Refs. 

For example, Seyfert galaxies possess a dusty torus of gas at distances intermediate 
between the BLR and NLR (narrow line region). An observer whose line of sight 
to the black hole intercepts this torus would see a heavily reddened (or completely 
extinguished) BLR and central continuum radiation but an unreddened NLR. This 
would be identified with a Seyfert 2 galaxy. If the line-of-sight does not intercept the 
torus, the central regions of the nucleus can be observed directly, leading to a Seyfert 
1 classification. Radio loud quasars are then objects in which the line-of-sight is close 
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to the jet cone of the source. In the cases in which we are not directly looking into 
the jet cone -blazars where relativistic effects produce highly variable and continuum 
dominated emission- emission from the BLR can be observed. Objects with larger 
inclinations have a less dominant central continuum flux, resulting in Fanaroff- Riley II 
(FRII) galaxies. If the torus surrounding the black hole obscures the BLR, a narrow 
line radio galaxy (NLRG) can be observed. It is not clear how FRI radio galaxies fit into 
such a scheme. Clearly, some (as yet unknown) physical mechanism, probably related 
to source power, produces different radio morphologies in FRI and FRII sources. 
Some bl azars may be beamed FRI objects, but there is a lack of broad-line FRI radio 
galaxies 214| . This make the classification within the unified scheme harder to achieve. 



2. Radiogalaxies 



FRII galaxies |219| are the largest known dissipative objects (non-thermal sources) in 
the Universe. Localized regions of intense synchrotron emission, known as "hot spots" , 
are observed within their lobes. These regions are presumably produced when the 
bulk kinetic energy of the jets ejected by a central active nucleus (supermassive black 
hole -|- accretion disk) is reconverted into relativistic particles and turbulent fields at 
a "working surface" in the head of the jets |22fll |. Specifically, the speed Vh with which 
the head of a jet advances into the intergalactic medium of particle density ng can be 
obtained by balancing the momentum flux in the jet against the momentum flux of 
the surrounding medium. Measured in the frame comoving with the advancing head, 
fh ~ [1 + (^e/'^j)^'^^]~\ where rij and fj are the particle density and the velocity 
of the jet flow, respectively, fj > fh for rie > nj, and the jet will decelerate. The 
result is the formation of a strong coUisionless shock, which is responsible for particle 



reacceleration and magnetic field amplification [221|. The acceleration of particles up 
to ultrarelativistic energies in t he ho t spots is the result of repeated scattering back 
and forth across the shock front |222| . Dimensional arguments suggest that the e nergv 
density per unit of wave number of MHD turbulence is of the Kolmogorov type 
and so for strong shocks the acceleration time for protons is ||224] 
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The FaranofF-Riley classification is based on one parameter, i?FR, the ratio of the distance between 
the regions of highest surface brightness on opposite sides of the central galaxy to the total extent 
of the source. Objects with R-pK < 0.5 are classified as FRI, whereas those with i?FR > 0.5 are 
classified as FRII. It is found that the brighter sources are all FRII class, although the distinction 
between classes is not clear cut in luminosities (for further details see page 220 of Ref. 218]). 
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where /3jet is the jet velocity in units of c, u is the ratio of turbulent to ambient magnetic 
energy density in the region of the shock (of radius R), and B is the total magnetic 
field strength. The acceleration process will be efficient as long as the energy losses 
by synchrotron radiation and photon-proton interactions do not become dominant. 
The subtleties surrounding the conve rsion of a particle kine t ic en ergv into radiation 
provide ample material for discussion IS H, H, E3, Q . The proton blazar 
model relates 7-ray emission to the development of electromagnetic cascades triggered 
by s econdarv photqmeson p roducts that cool instantaneously via synchrotron radia- 



tion 



iii l225[ liifil . I227I . l228 | . The synchrotron loss time for protons is given by jis 



T-syn ~ 



(Tt ml r i?2 



(31) 



where nip 



m. 



v 



(Tt and F are the electron mass, proton mass, Thomson cross section, 
and Lorentz factor, respectively. The characteristic single photon energy in synchrotron 
radiation emitted by an electron is 



1/2 



heE^B 



2 TT 



5.4 X IQ-'^B 



(32) 



For a proton this number is {mp/me)^ ~ 6 x 10^ times smaller. High energy 7-ray pro- 
duction through proton synchrotron radiation requires very large, (9(100 G), magnetic 
fields. Consi derin g an average cross section a-yp for the three dominant pion-producing 



interactions |231[, "fp — >■ pir^ ,'yp 



niT 



,7P 



pvT+vr , the time scal e of the energy 



losses, including synchrotron and photon interaction losses, reads 
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(33) 



where a stands for the ratio of photon to magnetic energy densities and A gives a 
measure of the relative strength of 7p interactions versus the synchrotron emission. 
Note that the second channel involves the creation of ultrarelativistic neutrons (but 
F„ < Fp) with mean free path in the observer rest frame given by An = FnCr„, where 



Tr,. ~ 



900 s, is the neutron lifetime. Since A„ > A„ for F„ < F 



p max; 



such neutrons can 

readily escape the system, thereby modifying the high end of the proton spectrum. 
Biermann and Strittmatter |222| have estimated that A ^ 200, almost independently 
of the source parameters. The most energetic protons injected in the intergalactic 
mediu m w ill have an energy that can be obtained by balancing the energy gains and 
losses |l25l | 

= 1.4 X 10^ B-t Aef u'^' Rlli' (1 + Aa)-''' , (34) 



where i?kpc = -R/lkpc. 
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For typical hot-spot conditions {B ~ 300 /iG, u ~ 0.5, and Pjet ~ 0.3) and assuming 
that the magnetic field of the hot spot is limited to the observable region, one obtains 
E < 5 X 10^° eV for a < 0.1 232l| .-^^ Particles can also attain ultrahigh energies 
{E > 10^° eV) within the jets or the AGNs themselves. For instance, the knot A 
in the M87 jet, with a length scale /gy ~ 2 x lO^'^ cm, has a magnetic field str ength 
~ 300 fiG [233J. Typical AGN sizes are /agn ~ 10^^ cm, and 5agn ~ 1 G 234 |. 
Observational evidence suggests that in the jets a <^ 1, whereas a ~ 1 for AGNs 222 1. 



3. Cen A: The source of most UHECRs observed at Earth'; 



Centaurus A (Cen A) is the nearest active galaxy, ~ 3.4 Mpc |235l |. It is a complex 
FRI radio-loud source identified at optical frequencies with the galaxy NGC 5128. 
Different multi-wavelength studies have revealed that it is comprised of a compact 
core, a jet also visible at X-ray frequencies, a weak counterjet, two inner lobes, a 
kpc-scale middle lobe, and two giant outer lobes. The jet would be responsible for the 
formation of the northern inner and middle lobes when interacting with the interstellar 
and intergalactic media, respectively. There appears to be a compact structure in the 
northern lobe, at the extrapolated end of the jet. This structure resembles the hot 
spots such as those existing at the extremities of FRII galaxies. However, at Cen A, 
it lies at the side of the lobe rather than at the most d istan t northern edge, and the 
brightness contrast (hot spot to lobe) is not as extreme |236| . 

Low resolution polarization m easur ements in the region of the suspected hot spot give 
magnetic fields as high as 25 /xG [236^ . However, in certain regions where measurements 
at both high and low resolution are available, the S-field amplitude at high resolution 
can be seen to be twice that at low reso lution. The higher resolution can reveal 



amplifi cation in the post-shock region |237| . yielding 5-fields possibly as high as 50 



60 /xG 238L 1239^ . T he ra dio- visible size of the hot spot can be directly measured from 
the large scale map 24ol |. giving i?HS — 2 kpc. The actual size can be larger by a factor 
~ 2 because of uncertainties in the angular projection of this region along the line of 
sight. Then, if the magnetic field of the hot spot is confined to the visible region, the 
limiting energy imposed by the Hillas' criterion is -Emax ~ 10^°'^ eV. 

Estimates of the radio spectral index of synchrotron emission in the hot spot and 



The shock structure in hot spots is hkely to be much more extended than the visibl e region in the 
non-thermal radioemission, as suggested by magnetohydrodynamical modeling |232l |. 
For example, an explanation of the apparent absence of a counterjet in Cen A via relativistic 
bearning suggests that the angle of the visible jet axis with respect to the line of sight is at most 
36° I233, which could lead to a doubling of the hot spot radius. It should be remarked that for a 
distance of 3.4 Mpc, the extent of the entire source has a reasonable size even with this small angle. 
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the observed degree of linear polarization in the same region suggests that the ra tio o f 
turbulent to ambient magnetic energy density in the region of the shock is m ~ 0.4 



241 



The i et ve locity is model dependent: possible values range from ~ 500 km s ^ to 



0.99 c |236l |. For FRI galaxies, the ratio of photon to magnetic energy densities, a, is 
expected to be -C 1. Now, by replacing these numbers into Eq. (jHll), one can easily see 
that Cen A can accelerate particles to energies > 10^" eV, with a maximum attainable 
energy set by the Hillas' criterion. 

Recent observations of the 7 ray flux for energies > 100 MeV by EGRET |243l | 
allow an estimate ~ 10^^ erg s~^ for the source. This value of is consistent 
with an earli er observation of photons in the TeV-range during a period of elevated 



activity 244| . and is considerably smaller than the estimated bolometric luminosity 



Lboi ~10 erg s |235l |. Data across the entire 7 ra y bandwidth of Cen A is given in 



R,ef . I245I 1 . reaching energies as high as 150 TeV 24fil |. though data at this energy await 



confirmation. For values of B in the /iG range, substantial proton synchrotron cooling 
is suppressed, allowing the production of high energy electrons through photomeson 



processes. The average energy of synchrotron photons scales as ^ 0.29-^^ |247 
With this in mind, it is straightforward to see that to account for TeV photons, Cen 
A should harbor a population of ultra-relativistic electrons with ii^ ~ 6 x 10^^ eV. We 
further note that this would require the presence of protons with energies between one 
and two orders of magnitude larger, since the electrons are produced as secondaries.^^ 



Th ere a re plausible physical arguments 228L 1248^ as well as some observational rea- 
sons 249l | to believe that when proton acceleration is being limited by energy losses, 
the CR luminosity Lcr ~ L^. Defining e, the efficiency of UHECR production com- 
pared to high energy 7 production - from the above, e ^ 1 - and using equal powe r 
per decade over the interval (-Emin, -^max), the source luminosity is found to be |250j 

(iiVg+" _ 6.3eL4i lO^W/s ^^^^ 

dEdt ln(£;max/^min) 

where L41 = luminosity of Cen A/lO^^erg s"^ and the subscript "0" refers to quantities 
at the source. 

For fiducial values, B = 0.5 ^G, £ = 0.5 Mpc, the diffusive distance traveled by CRs 
with E = 10^^ eV, is ctd = 50 Mpc ^ c? = 3.4 Mpc. Moreover, one can easily check 



It was recently suggested |242| | that the 25 kpc jet (with B ^ 10 /iG) of Cen A could be another 
promising region for proton acceleration to ultra high energies. 

Note that th e rec eived radiation is negligibly affected by interactions with the various radiation 
backgrounds j22d |. 

Consecutive factors of ~ 2 energy loss occur in the processes pj — > iV7r°, tt" — > 77, 7 e+e~. 
Ea. H32() then implies proton energies of ^ 10^° eV for 100 TeV photons. 
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that for 3.4 Mpc the diffusion time of any proton with energy above the photopion 
production threshold is always less than the GZK-time, and consequently energy losses 
can be safely neglected. This implies that the density of protons at the present time t of 
energy ii^ at a distance r from Cen A (which is assumed to be continuously emitting at 
a constant spectral rate (INq^^ /dE dt from time ton until the pres ent) can be obtained 



by solving the Kolmogorov-diffusive-equation, and is found to be 251 



rfn(r, t) dNT 1 r e-V^^ft-*') _ dNr 1 ... 



dE dEdt [A'kD{E)Y/^ Jt^^ {t-t'f/^ dEdtATxD{E)r 

where D[E) is the diffusion coefficient given in Eq. ((Tj), x = ADTo^/r"^ = Ton/r/), 
^on — t ton; and 

I{x) = ^r^e-- . (37) 

For Ton — * cxD, the density approaches its time- independent equilibrium value rieq, while 
for Ton = T"D, ri/riea, = 0.16. 

To estimate the power of Cen A, one can evaluate the energy-weighted approximately 
isotropic proton flux at 1.5 x 10^^ eV, which lies in the center of the fiat "low" energy 
region of the spectrum, 



In Eq. (j58jl we have used the fiducial values of B and ^ as given in the previous 
paragraph, and set -Emin = 1 x lO'^^ eV, -Emax = 4 x 10^° eV. As noted by Farrar 
and Piran |250| . by stretching the source parameters the "low" energy fiux from Cen 
A could be comparable to that of all other sources in the Universe. To this end, 
first fix eT4i I = 0.40, after comparing Eq. to the observed CR-fiux by AGASA: 
E3Jobs(^) = 10^^-^ eV^ m-2 s-i sr-i Next, 6^41 ~ 1, determines I ~ 0.40, and 
conseque ntly the r equired age of the source Ton to be about 400 Myr, which appears 
plausible |22l[E49l |. To maintain fiux at the "ankle" for the same Ton, one requires an 
approximate doubling of Lqr at 5 x 10^^ eV. Because of the larger diffusive time delay 
at this energy, this translates into an increased luminosity in the early phase of Cen 
A. From Eq. (j32j) . the associated synchrotron photons are emitted at energies < 30 
MeV. The increase in radiation luminosity in this re gion is no t inconsistent with the 
flattening of the spectrum observed at lower energies 

Having identified Cen A to plausibly be a powerful source of UHECRs, we now 
explore whether 5-field deflections provide correct directional properties, i.e., sufficient 
isotropy. This can be found by computing the incoming current flux density DVu as 
viewed by an observer on Earth, and one flnds for a continuously-emitting source a 
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distribution ~ (1 + acos6') about the direction of the source, where 6 is the angle to 
the zenith and 



a 



2D{E) I' 



X] 



cr 



—= [ du 

^ -ll/x 



u e 



(39) 



with X = Ton/Tc, and I as defined in Eq. ^7\i |25]| . For our choices of B and i, and 
Ton = 400 Myr, we find for E = 10^^ eV {E = 10^° eV) that a = 0.04 (a = 0.07). This is 
in complete agreement with the upper bounds on dipole anisotropies recently reported 
by HiRes Collaboration [t^ . One caveat is that the large defiection angle of the highest 
energy Fly's Eye even t with respect to the line of sight to Cen Amust be explained as 
a 2a fiuctuation 254l |. Additionally, Monte Carlo simulations 255l | show the predicted 



auto-correlation function is not consistent with the clustering at small scale reported 
by AGASA Collaboration 
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Therefore, if the hypothesis of CR pairing proposed by 
AGASA Collaboration is confirmed by future data, it will constitute a serious objection 
to the model outlined above. On the other hand, an interesting observational feature 
for a Cen A origin of UHECRs is the possible detection of neutrons, which at the hi ghes t 
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energies could survive decay and produce a spike in the direction of the source 
The estimated event rate at PAO is about 2 direct events per year, against negligible 
background. Thus, in a few years of running, the hypothesis of Cen A as the source of 
most UHECRs observed at Earth can be directly tested. 



4. M87: The end of all roads? 



M87 is a giant radi o ga laxv for which there has been a recent report of a TeV 
excess at a level of 4cr [256]. It is also expected to be a source for GLAST, having an 
EGRET upper limit of 2.8 x 10^^ photons cm^^ s~ ^ abo ve 100 MeV (Reimer, private 
communication, see *o the limit topo=ed in Ref. 0), and comparable theoretical 
fiux predictions 258, 259^ . 

M87 was thought as a ^high-energy CR emitter since quite long ago [2& 
a distance of 16.3 Mpc 
(/ = 282°, h = 74°) m 
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At 

262| . it is the dominant radio galaxy in the Virgo cluster 



The emission of synchrotron radiatio n with a s teep cutoff 



at frequencies about 3 x 10 Hz from its radiojets and hot spots [26J, l265j implies an 
initial turbulence injection scale having the Larmor radius of protons at 10^^ eV. 

The major difficulty with a M87 generation of UHECRs is the observation of the 
nearly isotropic distribution of the CR arrival directions. One can again argue that 
the orbits are bent. However, the bending cannot add substantially to the travel time, 
otherwise the energy would be GZK-degraded. An interesting explanation to overcome 
this difficulty relies on a Galactic wind, akin the solar win d, that would bend all the 
orbits of the highest energy CRs towards M87 \26(i l267l |. Indeed, it has long been 
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ex pecte d that such a kind of wind is active in our Galaxy |26iE6i,E23- In the analysis 
of |266| . it was assumed that the magnetic field in the Galactic wind has a dominant 
azimuthal component, with the same sign everywhere. This is because in a spherical 
wind the polar component of the magnetic field becomes negligible rather quickly, 
decaying like and thus the azimuthal part of themagnetic field quickly becomes 
dominant, with i?^ ~ sin 6'/r in polar coordinates 27ll |. Under these considerations 
one is left with two degrees of freedom: the strength of the azimuthal component at 
the location of the Sun, and the distance to which this wind extends . Re cent estimates 
suggest that the maguetie field strength near the Sun is ~ 7 H. The second 
parameter is more uncertain. Our Galaxy dominates its near environment well past our 
neighbor, M31, the Andromeda galaxy, and might well extend its sphere of influence 
to half way to M81. This implies an outer halo wind of ~ 1.5 Mpc. With this in mind, 
the mean flight time of the protons in the Galaxy is ~ 5.05 x 10^ yr ^ r^, the time 
for straight line propagation from M87 (Medina Tanco, private communication). The 
directions where the 13 highest energy CR events point towards when t hev l eave the 
halo wind of our Galaxy is consistent with an origin in the Virgo region |26q : (i) for 
CR protons, except for the two highest energy events, all other events can be traced 
back to within less than about 20° from Virgo; (ii) if one assumes that the two highest 
events are helium nuclei, all 13 events point within 20° of Virgo. Arguably, the super- 
Galactic plane sheet can focus UHECRs along the sheet. Hence, the particles would 
arrive at the boundary of our Galactic wind with the arri val d irections described by 



2721. This would allow a 



an elongated ellipse along the super-Galactic plane sheet 
bending of 20° to be accomodated. 

Additionally, in order to account fo r mo st of the CRs observed above the ankle, 
the power requirement of Virgo cluster H needs a fine-tuning of the source direction 



relative to the symmetry axis of the wind, so as to turn on magnetic lensing effects |274 
In such a case, M87 could be as high as > 10^ times more powerful than if unlensed at 
en ergie s below E/Z ~ 1.3 x 10^° eV. Criticisms of this model |275l | have been addressed 
in Em. 



5. Other powerful nearby radiogalaxies 

Apart from Cen A (which would provide the most energetic particles detectable on 
Earth), the CR-sky above PAO, if populated by radiogaxies, should be dominated by 
Pictor A (a str ong s ource with a flat ra dio s pectrum) which would contribute with the 
larger CR flux |232| . and PKS 1333-33 |277| . Other two southern candidates would be 
Fornax A {z = 0.057) and PKS 2152-69 {z = 0.027), which could provide contributions 
to t he CR fl ux above the cutoff. For other powerful sources and their properties 



see 
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There are two additional EGRET sources, one of them at high latitude, for which 
a possible radio galaxy counterpar t ha s been suggested. One such source is 3EG 
J1621+8203 (/ = 115.5°, 6 = 31.8°) H- 3EG J1621+8203 observations in individual 
viewing periods yielded near-threshold detections by EGRET, as for Gen A. However, 
in the cumulative exposure, it was clearly detected and the measured flux above 100 
MeV was 1.1 x 10^'' photon cm~^ s~^. The photon spectral index for this source is 
2.27±0.53, steeper than the usual blazar-like spectrum. Mukherjee et al. |279l | analyzed 
the X-ray and radio field coincid ent w ith 3EG J1621-(-8203 and concluded that NGG 
6251, a bright FRI radio galaxy 214| at a redshift of 0.0234 (implying a distance 91 
Mpc for Hq = 75 km s~^ M pc'-^ j , and the parent galaxy of a radio jet making an angle 
of 45° with the line of sight 280l |. is the most likely counterpart of the EGRET source. 
With this identification, the implied 7-ray luminosity is also a factor of 10~^ below that 
typical of blazars. Gompared with Gen A, the greater distance to NGG 6251 could, 
perhaps, be compensated by the smaller angle between the jet and the line of sight. 

Gombi et al. 28l| have also recently reported the discovery of a new radio 
galaxy, J1737— 15, within the location error box of the low- latitude 7-ray source 
3EG J1735— 1500, whose photon index is F = 3.24 ± 0.47. The radio galaxy mor- 
phology at 1.4 GHz is typical of the double-sided FRII. The integrated radio flux 
is 55.6 ± 1.5 mJy at 1.4 GHz, the source is non-thermal and it is not detected at 
4.8 GHz. Using the relation between approaching and receding jets: Sa.ppr/Sj.ec = 
{I + j3 cos^^/1 — (3 cos 6*) ~" , as well as the radio fluxes of each jet component, a view- 
ing angle in the range 79° — 86° for a velocity (3 = v/c between 0.3 and 0.9 and a = —1 
is derived. Depending on the jet and ambi ent m edium parameters, most double-sided 



radio sources have sizes below ~ 300 kpc j22l|. In the case of J1737— 15, and using 



standard Friedmann-Robertson- Walker models, this size translates into a possible dis- 
tance smaller than 350 Mpc. If 3EG J1735— 1500 is indeed the result of 7-ray emission 
in J1737— 15, the intrinsic luminosity a.t E > 100 MeV, at the distance quoted, should 
then be less than 2 x 10^"^ erg s~^, also several orders of magnitude smaller than that 
of blazars. If both radiogalaxies are closer than 100 Mpc, they could also be relevant 
acceleration sites of the observed UHEGRs. 



6. Correlations of UHECRs with QSOs, BL LACs, and EGRET sources 

Since an alignment beyond random expectations between UHEGRs and QSOs would 
certainly constitute a great discovery, the possible correlation between UHEGRs and 
QSOs was subject to a great deal of scrutiny. In the spring of 1998, Farrar and Bier- 
mann pointed out the existence of a directional correlation between compact radio- 
QSOs and UHEGRs: all events at the high end of the spectrum observed by that time, 
with energy at least la above 10^^'^ eV, were aligned with high redshifted quasars, a 
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phenomenon with a chance probabihty of occurrence less than 0.5% [57[. Since then, 
this correlation has been analyzed several times. Hoffman stated that one of the 5 
events used in the Farrar and Biermann's study, the highest energy event observed by 
the Fly's Eye experiment, should not be included in the UHECR sample under anal- 
ysis, because this very same event was considered to introduce the hypothesis |282 . 
Without this event, the positive alignment with random background probability is 
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increased to < 3%, in any case small enough as to be plausibly significant ^ ^ 

ing an updated event list (twice the size of th e pr evious) from the Haverah Park [i^ 
and the AGASA experiments, Sigl et al. |283l | showed that the statistical signifi- 
cance of th e ali gnment is lowered to 27%. Other authors, however, favored the earlier 
alignment |284| . but their correlation signal comes from events with large uncertainty 
both in energy and in position: they considered events from the SUGAR experiment, 
although it is not clear whether all these events are above the GZK cutoff. Not with- 
standing, after the Haverah Park energy estimate, have been re-a=.e.=ed the 
original correlation has to be dropped altogether: for the cosmic rays in question, the 
energy of the 2 events observed by this array with incident zenith angle < 45°, that 
was previously quoted as > 10^^'^ eV at la, is now shifted ~ 30% downwards, below 
the energy cut chosen by Farrar and Biermann. Hence, independently of the statistical 
test used, when considering only the highest energy (> 10^^'^ eV at lex) events the 
correlation between UHECRs and QSOs is consistent with a random distribution at 
the la level. 



Tinyakov and Tkachev |285L l28fiL l287| reported a correlation between the arrival 
directions of UHECRs and BL Lacs. Specifically, the (22) BL Lacs chosen were those 
identified as such in the (9th-Edition) Veron-Cetty and Veron (2000) |288l | catalogue 
of Quasars and Active Galactic Nuclei, with redshift 2; > 0.1 or unknown, magnitude 
m < 18, and radio flux at 6 GHz Fq > 0.17 Jy. This analysis proposed no energy buffer 
against contamination by mis-measured protons piled up at the GZK energy limit. 
The evidence supporting their claim is based on 6 events reported by the AGASA 
Collaboration (all with average energy < 10^^'^ eV), and 2 events recorded with the 
Yakutsk experiment (both with average energy < 10^^'^ eV), which were found to be 
within 2.5° of 5 BL Lacs contained in the restricted sample of 22 sources. The chance 
probability for this coincidence set-up was claim to be 2 x 10~^. Here also the data set 
used to make the initial assertion is also being used in the hypothesis testing phase. 
What is further subject to critique, is that the imposed cuts on the BL Lac catalogue 
were chosen so as to maximize the signal-to-noise ratio, compensating a posteriori the 



The CR sample of Tinyako v and Tkachev consists of 26 events measured by the Yakutsk experiment 
with energy > lO^^ -^^ eV |289j . and 39 events measured by the AGASA experiment with energy 



> 1019-68 eV 
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different cut adjustments by inclusion of a penalty factor j290|. Without such arbitrary 
cuts, the significance of the correlation signal is reduced at the la level. Not to anyone's 
surprise, even in acc eptance of this approach, the estimated value of the penalty factor 
is subject to debate 287, 290^ . 

Recently, i n ord er to test the hypothetical correlation between UHE CRs and BL Lacs, 



Torres et al. 
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performed a blind analysis using the Haverah Park ^29^ and Volcano 
Ranch 2931 ] data samples. Such an analysis shows no positional coincidences between 
these two samples up to an angular bin > 5°, an angular scale that is well beyond 
the error in arrival determination of these experiments (~ 3°) 'tt']. On the basis of 
the strongly correlated sample analyzed by Tinyakov and Tkachev, one expects the 
distribution describing the correlation between the set of BL Lacs and any UHECR 
data-set with 33 entries to be Poisson with mean ^ 4.06. This implies a 2a deviation 
effect. Alternatively, the 95% CL interval of the distribution which samples the corre- 
lation between the BL Lacs and CRs recorded by Volcano Ranch + Haverah Park is 
(0, 3.09) (294| . so that the probability tomeasure the expected mean value ~ 4.06 is 
^ 5%. With this in mind, Torres et al. |29lj conclude that the 8 coincidences found in 
the Tinyakov and Tkachev's analysis do not represent a statistically significant effect. 

Additionally, Gorbunov et al. |295| claimed that a set of 7-ray loud BL Lacs can 
be selected by intersecting the EGRET, the UHECR, and the BL Lac catalogs (all 
conveniently cut). The only requirement Gorbunov et al. considered for an object 
(here, a BL Lac) to be physically associated with an EGRET source is that the angular 
distance between the best estimated position of the pair does not exceed 2 x Rg^, where 
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i?95 is the 95% confidence level contour of the EGRET detection. Torres et al. 
pointed out that identifying EGRET sources with BL Lacs (or any other object) just by 
positional pairing within twice the EGRET error grossly underestimates the goodness 
of existing 7-ray data. At this stage, it is worth recalling the reader that the typical 
i?95 radius for EGRET sources is 0.5-1°. One can then argue that if the confidence 
contours have any significance at all, a source should appear beyond the 95% contour 
only a f ew p ercent of the time. Working with 114 EGRET sources above |6| > 10°, 
Punsly (296l | have estimated the number of random coincidences as a function of the 
field radius: ~ 2 (10) quasars with more than 1 Jy of 5 GHz flux are expected to 
correlate by random chance if the size of the typical EGRET angular uncertainty is 
0.7° (1.7°), see Fig. El 

In our opinion, available statistics on the arrival directions of the UHECRs reveals 
no significant correlations above random with BL Lacs nor with any other type of 
quasars, including EGRET blazars. 
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FIG. 6: The expected distribution of radio-loud quasars (louder than 0.5 Jy at 5 GHz) to 
occur by random chance as a function of the distance from the center of the field for a sample 
of 114 EGRET detections. Points represent the number of 7-ray detections for which the 
counterparts are beyond the 95% confidence contour. The dotted curve are the boundaries of 
the 68% confidence band for the hypothesis that the radio sources are randomly distributed 
in the EGRET detection fields. Adapted from Punsly (1997). The number of sources whose 
possible counterpart are beyond the 95% confidence contour is compatible with the chance 
expectation. 



E. Remnants of quasars 

1. What is a quasar remnant and how would they accelerate particles? 

Interestingly, the absence of powerful radio emitting objects in the direction of sev- 
eral UHECRs led some colleagues to think that dead, faint objects, yet ones sufficiently 
active as to accelerate particles up to relativistic energies, are responsible for the UHE- 
CRs observed. S uch is the idea behind the concept of quasar remnants (QRs) as 
UHECR emitters 2971 l298l |: a spinning supermassive black hole, threaded by mag- 
netic fields generated by currents fiowing in a disc or torus, induces an emf which, if 
vaccum breakdown is prevented and in the absence of severe energy losses, accelerate 
a particle near the full voltage. 

Although in the present epoch there is a paucity of luminous quasars (L > 10'^'' ergs 
s~^), those which appear at large redshifts, the expected local number of dead quasar s 
associated with the same parent population is expected to be large 



299, M, 301 
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Supermassive black holes are now to be found in the relatively dormant nuclei of giant 
elliptical galaxies, generally regarded as QRs.^^ The compact dynamo model has b een 
proposed as a n atural mechanism for accelerating CRs in such environments |297 |. In 
this model jsO^, if B is the ordered poloidal field near the hole, V ~ aB, where a is 
the hole's specific angular momentum (e.g., a = M for an extreme Kerr hole of mass 
M). In an appropriate astrophysical scaling |304l | 

V ^9x 1020(a/M)54M9 V, (40) 

where B^ = 5/(10^G) and Mg = M/(1O^M0). Assuming that the energy density of 
the magnetic field near the e vent horizon is in equipartition with the rest mass energy 
density of accreting matter |305| . it is possible to introduce the accretion rate in the 

voltage expression; In terms of an Advection Dominated Accretion Flow (ADAF) 

model (e.g., Ref. jsO^), for example, 

^4 ~ 1.4M9-^M^/^ (41) 

where M is the accretion rate dM/dt in Mq yr"-*^, and then, the maximum emf {a/M ~ 
1) is 

V r^l.2x lO^^M^/^ V. (42) 

This potential is not, however, the maximum obtainable energy. 

The rate of energy loss through curvature radiation by a particle of energy E = mc^T 
can then be expressed as 



3 p2 



^=0-^' (43) 



see 



3071 ] for a detailed explanation. Here, p is the average curvature radius of an 
accelerating ion, assumed to be independent of the ion energy. The energy change per 
unit length of an accelerating ion having charge Z and mass rrii = prrip is given by 

de/ds = eZAV/h - P/c, (44) 

where h is the gap height. After integration from s = Oto h, the maximum acceleration 
energy is obtained 

E^,^ = 3 X 10'' pZ'/^Ml/'Bl/\p^h/RlY/^ eV. (45) 
Consequently, only a fraction 

V = 0.1pM,'/\ZB,r'/\p/Rgf\h/R,r'/'; < 1, (46) 



Indeed, the term "quasar remnants" was introduced by Chokshi and Turner (302| to describe the 
present-epoch population of dead quasars harboring supermassive black hole nuclei. 
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of the potential energy available will be released as UHECRs; the rest will be radiated 
in the form of curvature photons. 

For a proton, the suppression ratio is E/[e{V)\ ^ [{50Mg)-^/^B4'^/y^/\ where r 
deno tes the magnetic field curvature in units of the Schwarzschild radius {h ~ Rg) 
j298l |. For r ^ 1 and M ^ (0.1 — 10) M© yr~^, and using the previous equations, 

^max = (1.0 - 1.8) X lO^^Mgl/^ cV. (47) 

Heavier nuclei would reach higher energies, but are subject to photo-disintegration. 
For highly energetic protons, energy losses due to photo-pion production in collisions 
with ambient photons also becomes a relatively important effect. A lower limit to the 
radiation length (Amin) for the proton energy loss associated with photo-pion produc- 
tion is estimated by considering the population of target photons within the source 
region [i?(source radius) > 2GM/c^] at radio frequencies u > 360(r/10^^)^^ GHz is 
given by |298l | 

A/R = CTiR/{{Ka.p^)Q) (48) 
= (278/(ir(Tp^,^b))(i?/i?5)M9(Q/10^3 

where Cp^^^y, = crp^/10~^° cm^, Rs is the Schwarzschild radius (twice the gravitat ional 
radius Rg = GM), K = (i? (loss) (initial) is the inelasticity in a single collision |308l | 
and Q is the core emission rate (photons s~^) for electromagnetic radiation aX v > 360 
GHz, given by Q = J v~^Lydv where h is the Planck constant and Ly = AttD'^F^ 
or a source of spectral density F,^ at distance D. A useful approximation is that 
298l | [Kap^^^]^) = {Kap^i^^\,{.dQ / di')dp\ /Q < 120 /ib. Those quasar remnants with 
Amin > R are expected to successfully accelerate protons up to ~ -Emax- Noteworthy, 
the observed flux of CRs would apparently drain only a negligible amount of energy 
from the black hole dynamo, since replenishing the particles ejected at high energies 
(> lO^'' eV) would only require a minimal mass input: a CR luminosity of 10^^ ergs/s 
in such particles (if protons) corresponds to a rest mass loss rate of less than 1O~^M0 
in a Hubble time. 



2. Correlations of UHECRs with QRs 

The first analysis of possible correlations between UHECRs and QRs was carried 



out in Ref. 309|, where, imposing very restrictive selection criteria -those necessary 
for obtaining candidate objects providing the most favorable setting for a black hole 
based compact dynamo model of UHECR production- a gro up of galaxies from the 
Nearby Optical Galaxy (NOG) catalog of Giuricin et al. l310t| was a priori selected as 
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individual plausible sources of CRs.^^ It was found that nearby QR candidates present 
an above-random positional correlation with the sample of UHECRs. Surprisingly, 
this correlation appears on closer angular scales than that expected when taking into 
account the deflection caused by typically assumed intergalactic or Galactic magnetic 
fields. 26 

As one can see using Eq. (jHJ, scattering in large scale magnetic irregularities O (nG) 
are enough to bend the orbi ts of super-GZK protons by about 4 deg in a 50 Mpc 
traversal (see, how ever, Ref. |31lj ). The CR angular offsets observed for the quasar 
remnants in Ref. 309l | are much smaller than 6 whereas, for a variety of assumed 
magnetic field scenarios, 6 is often substantially larger than the estimated AGASA 
measurement error of 1.6 degrees. Should the apparent clustering of correlated pairs 
be supported by future data, perhaps a viable scenario under which this could occur 
is that in which the intergalactic medium between Earth and the three apparently 
'clustered' quasar remnants is sufficiently different from the intergalactic medium in 
front of the remaining nine objects that are much more uniformly distributed on the 
accessible sky. This appears at least plausible judging from the 100 micron map of the 
region shown in Fig. [7| 

For the deflection of an energetic f 60 EeV) proton in traversing 34 Mpc (the mean 
of the QRs distances in Table 3 of [309|) to be less than a degree, the magnetic fleld 
must fulflU B <2x 10"^° ^-^/^ q ^j^j^^ appears to be not drastically different from 
the canonical nG B fleld and the usually assumed coherence length ~ 1 Mpc. Also, in 
some directions, the magnetic fleld of our own galaxy could well lead to a deflection 
of up to several degr ees for primaries with energies below 60 EeV; see, for instance. 
Table 1 of Ref. (l33| . A recent study of this issue was presented by Alvarez- Muniz 



et al. 312] • However, a possible fllamentary topology of the Galaxy's magnetic fleld 
would likely allow some directional windows, albeit narrow, where the deflection of an 

UHECR could be much less than typical. 

In a very recent paper, Isola et al. 1313^ have studied predictions for large and 
small scale UHECR arrival direction anisotropies in a scenario where the particles 
are injected with a mono-energetic spectrum (all particles coming from a source are 
emitted with the maximal energy of acceleration for that source as derived in that 
same paper) by a distribution of QRs. They flnd that the sample of (37) QRs they 
considered is distributed too anisotropically to explain the isotropic ultra high energy 
CR flux except in the case where extragalactic magnetic flelds of ~ O.l/iG extend 



The latter catalog is a complete magnitude- limited (corrected blue total magnitude B < 14), 
distance-limited (redshift z < 0.02) sample of several thousand galaxies of latitude |6| > 20°. 
Deflections due to the magnetic fields would of course be avoided if the primary were a photon, 
generated in the neighborhood of the QR via an accelerated charged particle interaction. 
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FIG. 7: IRAS 100 micron view of the North Galactic Pole (white=high flux, dark blue=low 
flux), towards the region of several QRs that appear coinciding with UHE CRs. Small circles 
mark the 12 galaxies (candidate quasar remnants) in the sample of Ref. |309l |. Red circles 
mark galaxies perhaps associated with UHECRs; white circles mark galaxies without asso- 
ciated UHECRs. Note that the palusible UHECR sources tend to be found in directions 
of lower 100 micron flux, [UHECR-coinciding source directions present IRAS fluxes 0.614 ± 
0.022 MJy/sr whereas the directions towards those galaxies non-coinciding with UHECRs 
presnet 1.68 ± 0.63 MJy/sr.] To the extent that the 100 micron flux traces the Galactic dust 
and magnetic fleld, then CRs are likely to be better aligned with their sources in directions 
of low flux. Figure courtesy of Tim Hamilton. 

over many Mpc. As statistical quantities for this analysis spherical multi-poles and 
the autocorrelation function were used. For a weak magnetic field, of order InG, the 
predictions appear to be inconsistent with the observed distribution of arrival directions 
of UHECRs, because the magnetic field is too weak to isotropize the distribution coming 
from a lim ited number of non-uniformly distributed sources, as already pointed out in 
Ref. [255^ . Isola et al. also found that the contribution from the farthest sources is 
completely negligible even for this weak magnetic field. None of the objects that were 
found close to UHECR positions in the Torres et al. previous set were included in the 
Isola et al. sample. Little further progress regarding possible correlation of sources 
can be made until a much more larger set of UHECRs is recorded, something that will 
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have to wait to the operation of PAO. 



3. TeV emission from QRs 

A con comitant effect of UHECR emission from QRs is that, as shown by Levin- 
son j314l | , the dominant fraction of the rotational energy extracted from the black hole 
is radiated in the TeV band. He showed that the spectrum produced by the curvature 
radiation of a single ion will peak at an energy 

^^^ax = l.5T'hc/p (49) 

= 1.6 X 10-^E^ax/i-'(^54)^/'(V^9)'^' (50) 

= 5M^/\ZB^Y/\p^hyRlY/^TeV, (51) 

and is a power law I{E^) oc (E^Y^'^ below the cutoff. The overall spectrum of curvature 
photons would depend on the energy distribution of the accelerating particles, and is 
expected to be somewhat softer below the peak. For E^ax = 3 x 10^° eV and h Rg, 
^7 max ^ 50fi~\ZB4y/^ TcV. Then, provided that vaccum breakdown does not occur, 
and that TeV photons can escape the system, QRs should emit 7-rays. 

Indeed, it was recently noted by Neronov et al. |315| that the concomitant TeV 
radiation would be at a level sufficiently high as to be (for many combination of the 
system parameters) ruled out by bounds imposed by HEGRA/AIROBICC. 

Recent numerical simulations additionally suggest that the accretion process and 
magnetic field structure in the vicini tv o f the horizon can be non-stationary, owing 
to rapid magnetic field reconnection |31fil |. This would probably lead to appreciable 
complications of the model, as the location of the gap, the injection of seed particles 
into the gap, and, perhaps, the voltage drop across it might change with time. How 
should this affect the picture described above is unclear at present. 

F. Starbursts 

1. What are they? 

Starbursts are galaxies (sometimes, the term also refers only to particular regions of 
galaxies) undergoing a large-scale star formation episode. They feature strong infrared 
emission originating in the high levels of interstellar extinction, strong HH-region-type 
emission-line spectrum (due to a large number of O and B-type stars), and considerable 
radio emission produced by recent SNRs. Typically, starburst regions are located close 
to the galactic center, in the central kiloparsec. This region alone can be orders of 
magnitude brighter than the center of normal spiral galaxies. From such an active 
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region, a galactic-scale superwind is driven by the collective effect of supernovae and 
particular massive star winds. The enhanced supernova explosion rate creates a cavity 
of hot gas (~ 10^ K) whose cooling time is much greater than the expansion time scale. 
Since the wind is sufficiently powerful, it can blow out the interstellar medium of the 
galaxy, preventing it from remaining trapped as a hot bubble. As the cavity expands, 
a strong shock front is formed on the contact surface with the cool interstellar medium. 
The shock velocity can reach several thousands of kilometers per second and ions like 
iron nuclei can be efficiently accelerated in this scenario, up to ultrahigh energies, by 
Fermi's mechanism |317| . If the supe r- GZ K particles are heavy nuclei from outside our 
Galaxy, then the nearby (~ 3 Mpc ^1^) starburst galaxies M82 {I = 141°, 6 = 41°) 
and NGC 253 (/ = 89°, b = —88°) are prime candidates for their origin. 



2. M82 and NGC253 



M82 is probably the best studied starburst galaxy, located at only 3.2 Mpc. The 
total star formation rate in the central parts is at least ~ 10 yr^^ [319]. The far 
infra red luminosity of the inner region within 300 pc of the nucleus is ~ 4 x 10^*^ 
|32o| . There are_~ Ix 10'^ Mq of ionized gas and ~ 2 x 10^ Mq of neutral gas in 
the IR so urce |32ol . l32ll |. The total dynamical mass in this region is ~ fl — 2) x 10^ 

Mr: ' ' 
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The main observational features of the starburst can be modelled with a 
Salpeter IMF extending from 0.1 to 100 Mq. The age of the starburst is estimated in 
~ (1 - 3) X 10^ yr ^2^. Around ~ 2.5 x 10^ Mq (i.e. ~ 36 % of the dynamical mass) 
is in the form of new stars in the burst 
with large numbers of early-type stars. 



32l|. The central region, then, can be packed 



322 



323 



3241). A 



NGC 253 has been extensively studied from radio to 7-rays (e.g. 
TeV detection was reported by CANGAROO 1325^ . but has been yet unconffimed by 
other experiments. More than 60 individual compact radio sources have been detected 
within the central 200 pc |326| . most of which are supernova remnants (SNRs) of only a 
few hundred years old. The supernova rate is estimated to be as high as 0.2 — 0.3 yr~^, 
comparable to the massive star formation rate, ~ O.lMe yr- 0^3. The central 
region of this starburst is packed with massive stars. Four young globular clusters near 
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the center of NGC 253 can account for a mass well in excess of 1.5x10^Mq 
Assuming that the star formation rate has been continuous in the central region for 
the last 10^ yrs, and a Salpeter IMF for 0.08-100 Mq, the bolometric luminosity of 
NGC 253 is consistent with 1.5 xlO^M© of young stars |328l |. Based on this evidence, 
it appears likely that there are at least tens of millions of young stars in the central 
region of the starburst. These stars can also contribute to the 7-ray luminosity at high 



energies 



1^ 



Physical, morphological, and kinematic evidence for the existence 
of a galactic superwind has been found for NGC 253 |33lj . Shock interactions with 
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low and high density clouds can produce X-ray continuum and optical line emission, 
respectively, both of which have been directly observed. 

A region about 1 kpc of the M82 galactic center appears to be a fossil starburst, 
presenting a main sequence stellar cutoff corresponding to an age of 100-200 Myr and a 
current average extinction of 0.6 mag (compare with the extinction of the central and 
current starburst region, 2.2 m ag) w hereas, nearby globular glusters age estimations 
are between 2 x 10* and 10^ yr [332]. It appears possible for this galaxy, then, that a 
starburst (known as M82 "B" ) of similar amplitude than the current one was active in 
the past. 



3. Two-step acceleration-process in starbursts 

The acceleration of particles in starburst galaxies is thought to be a two-stage process 



317[. First, ions are thought to be diffusively accelerated at single SNRs within the 
nuclear region of t he galaxy. Energies up to ~ 10^^~^^ eV can be achieved in this 
step (see, e.g. 1333^ ). Due to the nature of the central region, and the presence of the 
superwind, the escape of the iron nuclei from the central region of the galaxy is expected 
to be dominated by convection.^'' Collective plasma motions of several thousands of 
km per second and the coupling of the magnetic field to the hot plasma forces the CR 
gas to stream along from the starburst region. Most of the nuclei then escape through 
the disk in opposite directions along the symmetry axis of the system, being the total 
path travelled substantially shorter than the mean free path. 

Once the nuclei escape from the central region of the galaxy they are injected into 
the galactic-scale wind and experience further acceleration at it s terminal shock. CR 
acceleration at superwind shocks was first proposed in Ref. |335| in the context of our 
own Gala xv. T he scale length of this second shock is of the order of several tens of kpc 
(see Ref. j318|), so it can be considered as locally planar for calculations. The shock 
velocity Vsh can be estimated from the empirically determined superwind kinetic energy 
flux i?sw and the mass flux M generated by the starburst through: Eg^ = l/2Mv'^^. 
The shock radius can be approximated by r Vsh^, where r is the starburst age. 



The relative importance of convection and diffusion in the escape of the CRs from a region of disk 
scale height h is given by the dimensionless parameter, q = Vq h/n p, wh ere Vq is the convection 
velocity and kq is the CR diffusion coefficient inside the starburst 334]. When g ^ 1, the CR 
outflow is difussion dominated, whereas when g > f it is convection dominated. For the central 
region of NGC 253 a convection velocity of the order of the expanding SNR shells ~ 10000 km s~^, 
a sca le height /i ~ 35 pc, and a reasonable value for the diffusion coefficient Kf) ^ 5 x 10^^ cm^ s~^ 
|l93l |. lead to q ^ 216. Thus, convection dominates the escape of the particles. The residence time 
of the iron nuclei in the starburst results t^ES ^ h/Vo ~ 1 x 10^"'^ s. 



47 



Since the age is about a few tens of million years, the maximum energy attainable 
in this configuration is constrained by the limited acceleration time arising from the 
finite shock's lifetime. For this second step in the acceleration process, the photon 
field energy density drops to values of the order of the cosmic background radiation 
(we are now far from the starburst region), and consequently, iron nuclei are safe from 
photodissociation while energy increases to ~ 10^° eV. 

To estimate the maximum energy that can be reached by the nuclei, consider the 
superwind terminal shock propagating in a homogeneous medium with an average mag- 
netic field B. If we work in the frame where the shock is at rest, the upstream flow 
velocity will be vi (|vi| = Vsh) and the downstream velocity, V2. The magnetic field 
turbulence is assumed to lead to isotropization and consequent diffusion of energetic 
particles which then propagate acco rding to the standard transport theory [336| . The 



acceleration time scale is then 18871 1: tacc = ^ where k, is the upstream diffusion coef- 
ficient which can be written in terms of perpendicular and parallel components to the 
magnetic field, and the angle 6 between the (upstream) magnetic field and the direction 
of the shock propagation: k = k\\ cos^ 9 + kj_ sin^ 9. Since strong turbulence is expected 
from the shock we can take the Bohm limit for the upstream diffusion coefficient paral- 
lel to the field, i.e. k\\ = ^E/ZeBi, where Bi is the strength of the pre-shock magnetic 
field and E is the energy of the Z-ion. For the k± component we shall assume, fol- 
lowing Biermann 1888^ . that the mean free path perpendicular to the magnetic field is 
independent of the energy and has the scale of the thickness of the shocked layer (r/8). 
Then, = 1/8 r{vi — V2) or, in the strong shock limit, k_l = rv'l/12. The upstream 
time scale is tacc ~ r/(8t;i), r/^ivi = A/vf {E / (SZeBi) cos"^ 9 + rvf/Usin^ 9) . Thus, 
using r = ViT and transforming to the observer's frame one obtains 

E^a. ^ \zeBvlT ^ \zeB^T. (52) 
The predicted kinetic energy and mass fluxes of the starburst of NGC 258 derived 



from the measured IR luminosity are 2x 10^^ erg s~^ and 1.2 yr~^, respectively 818 
The starburst age is estimated from numerical models that use theoretical evolutionary 
tracks for individual stars and make sums over the entire stellar popu latio n at each 
time in order to produce the galaxy luminosity as a function of time |82Cl| . Fitting 
the observational data these models provide a range of suitable ages for the starburst 
phase that , in the case of NGC 258, goes from 5 x 10"^ to 1.6 x 10*^ yr (also valid for 



M82) j320|. These models must assume a given initial mass function (IMF), which 
usually is taken to be a power-law with a variety of slopes. Recent studies has s hown 
that the same IMF can account for the properties of both NGC 258 and M82 889l . 



Finally, the radio and 7-ray emission from NGC 258 are well matched by models with 
B ~ 50/x G [,828,] . With these figures, already assuming a conservative age r = 50 Myr, 
one obtains a maximum energy for iron nuclei of -E^ax > 3.4 x 10^° eV. 
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4- The starhurst hypothesis: UHECR-luminosity and correlations 



For an extragalactic, smooth, magnetic field of ~ 15 — 20 nG, diffusive propagation 
of nu clei below 10^'' eV evolves to nearly complete isotropy in the CR arrival direc- 



34nl . l34ll | Thus, we could use the rates at which starbursts inject mass, metals 



tions 

and energy into superwinds to get an estimate of the CR-injection spectra. Gener- 
alizing the procedure discussed in Sec. HID 31 -using equal power per decade over the 
interval 10^^'^ eV < E < 10^°'^ eV - we obtain a source CR-luminosity 



dE dt 



3.5elO^W/s (53) 



where e is the efficiency of ultra high energy CR production by the superwind kinetic 
energy flux. With this in rnind, the energy-weighted, approximately isotropic nucleus 
flux at 10^9 eV is given by H 

^'^(^) - (J,kE) 1Z° - " ' 10'»eI.eV^.-s-.-. ,54) 

where h = Im82 + -^ngc 253- To estimate the diffusion coefficient we used B^g = 15, 
•^Mpc = 0.5, and an average Z = 20. We fix 

eJ^ = 0.013, (55) 

after comparing Eq. ()54|) to the observed CR-flux. Note that the contribution of 
/m82 and /ngc 253 to critically depends on the age of the starburst. The relation 
"starburst-age/superwind-efficiency" derived from Eq. (j^ . leads to e ~ 10%, if both 
M82 and NGC 253 were active for 115 Mvr . The power requirements may be reduced 
assuming contributions from M82 "B" 340|. 

Above > 10^°'^ eV iron nuclei do not propagate diffusively. Moreover, the CR- 
energies get attenuated by photodisintegration on the CMB and the intergalactic in- 
frared background photons. However, the energy- weighted flux beyond the GZK-energy 
due to a single M82 flare 

E^J{E) = _^:^o^e-«*/56 ^ 2.7 x lO^^E^oee-^^/^^ eV^ m-^ s-^r-^, (56) 
[4:11 d)'^ dEo dt 



is easily consistent with o bser vation |340| . Here, R is the effective nucleon loss rate of 
the nucleus on the CBM |l07| . 

In the non-diffusive regime (i.e., 10^°'^ eV ^ E < 10^°'^ eV), the accumulated 
deflection angle from the direction of the source in the extragalactic -B-field is roughly 
10° ~ ^ ^ 20° |34ll |. The nuclei suffer additional deflection in the Galactic magnetic 
field. In particular, if the Galactic field is of the ASS type, the arrival direction of 
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FIG. 8: Left: Directions in Galactic coordinates of the four highest energy CRs at the 
boundary of the Galactic halo. The diamonds represent the observed incoming directions. 
The circles and arrows show the directions of neon and iron nuclei, respectively, before de- 
flection by the Galactic magnetic field. The solid line is the locus of incoming directions at 
the halo for other species with intermediate atomic number. The stars denote the positions 
of M82 and NGC253. The dashed lines are projections in the {l,b) coordinates of angular 
directions within 20° and 30° of the starbursts. Right: Curves of constant probabilities in the 
two-dimensional parameter space defined by the size of the c one a nd the minimum number 
of events originating within the resulting effective solid angle [342i l . 



the 4 highest energy CRs can be traced backwards to one of the starbursts 842l . 
Figure IHl shows the extent to which the observed arrival directions of the highest energy 
CRs deviate from their incoming directions at the Galactic halo because of bending 
in the magnetic field given in Eq. (|T!?|l . The incoming CR trajectories are traced 
backwards up to distances of 20 kpc away from the Galactic center, where the effect of 
the magnetic field is negligible. The diamond at the head of each solid line denotes the 
observed arrival direction, and the points along these lines indicate the direction from 
which different nuclear species (with increasing mass) entered the Galactic halo. In 
particular, the tip of the arrows correspond to incoming directions at the halo for iron 
nuclei, whereas the circles correspond to nuclei of neon. Regions within the dashed 
lines comprise directions lying within 20° and 30° degrees of the starbursts. It is seen 
that trajectories for CR nuclei with Z > 10 can be further traced back to one of the 
starbursts, within the uncertainty of the extragalactic deviation. 

The effects of the BSS configuration are completely different. Because of the aver- 
aging over the frequent field reversals, the resulting deviations of the CR trajectories 
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are markedly smaller, and in the wrong direction for correlation of current data with 
the starburst sources. We note that the energy-ordered 2D correlation distribution of 
the AGASA data is in disagreement with expectations for positively charged particles 
and the BSS configuration 13121]. 

We now attempt to assess to what extent these correlations are consistent with 
chance coincidence. We arrive at the effective angular size of the source in a two- 
step process. Before correcting for bias due to the coherent structure of the Galactic 
magnetic field, the deflections in the extragalactic and Galactic fields (regular and 
random components) may be assumed to add in quadrature, so that the angular sizes 
of the two sources are initially taken as cones with opening half-angles between 40° and 
60°, which for the purpose of our numerical estimate we approximate to 50°. However, 
the global structure of the field will introduce a strong bias in the CR trajectories, 
substantially diminishing the effective solid angle. The combined deflections in the / 
and b coordinates mentioned above concentrate the effective angular size of the source 
to a considerably smaller solid angle. As a conservative estimate, we retain 25% of this 
cone as the effective source size. A clear prediction of this consideration is then that 
the incoming flux shows a strong dipole anisotropy in the harmonic decomposition. 

Now, by randomly generating four CR positions in the portion of the sky accessible 
to the existing experiments (declination range 6 > —10°), an expected number of 
random coincidences can be obtained. The term "coincidence" is herein used to label 
a synthetic CR whose position in the sky lies within an effective solid angle QcR of 
either starburst. Q^s is characterized by a cone with opening half-angle reduced from 
50° to 24° to account for the 75% reduction in effective source size due to the magnetic 
biasing discussed above. Cosmic ray positional errors were considered as circles of 1.6° 
radius for AGASA. For the other experiments the asymmetric directional uncertainty 
was represented by a circle with radius equal to the average experimental error. The 
random prediction for the mean number of coincidences is 0.81 ± 0.01. The Poisson 
probability for the real result to be no more than the tail of the random distribution 
is 1%. Alternatively, we may analyze this in terms of confidence intervals. For the 4 
observed events, with zero background, the Poisson signal mean 99% confidence interval 



is 0.82 — 12.23 j294| . Thus our observed mean for random events, 0.81 ± 0.01, falls 
at the lower edge of this interval, yielding a 1% probability for a chance occurrence. 
Of course, this is not compelling enough to definitively rule out chance probability as 
generating the correlation of the observed events with the candidate sources, but it is 
suggestive enough to deserve serious attention in analyses of future data. 

Assuming an extrapolation of AGASA flux {E^ Johs{.E)) up to 10^°'^ eV, the event 



Because of constraints inherent in partitioning events among clusters, the distributions are very 
close to, but not precisely Poisson 76]. 
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rate at Pampa Amarilla^^ is given by 



dN 



A 



El ^ 



A 



1 

Ef 



E"^ 



5.3 yr" 



(57) 



where Ei = 10^°'^ eV and E2 = 10^^'^ eV. Considering a 5-year sample of 25 events 
and that for this energy range the aperture of PAO is mostly receptive to cosmic rays 
from NGC253, we allow for different possibilities of the effective reduction of the cone 
size because of the Galactic magnetic field biasing previously discussed. In Fig. |H1 
we plot contours of constant probabilities [P = 10~^, 10~^) in the two-dimensional 
parameter space of the size of the cone (as a fraction of the full 50° circle) and the 
minimum number of events originating within the resulting effective solid angle. The 
model predicts that after 5 years of operation, all of the highest energy events would be 
observed in the aperture described above. Even if 7 or 8 are observed, this is sufficient 
to rule out a random fluctuation at the 10~^ level. Thus, a clean test of the starburst 
hypothesis can be achieved at a very small cost: < 10~^ out of a total 10"'^ PAO 
probability budget |343l |. 



G. Luminous Infrared Galaxies 

1. Definition 

In carrying to the extreme the concept of a starburst, we find the powerful lumi- 
nous infrared galaxies (LIGs). At luminosities above 10^^ L©, LIGs (LpiR > 10^^ L©) 
are the dominant extragalactic objects in the local universe {z < 0.3) with such high 
lumi nosities. Some, having LpiR > 10^^ Lq, are the most luminous local objects (see 
I344I I for a revi e w) . These gal axies possess very large amounts of molecular material 



(e.g., |345L I346L I347L I348L l349| ). They present large CO luminosities, but also a high 
value for the ratio Lfir/Lco, both being about one order of magnitude greater than 
for spirals, implying a higher star formation rate per solar mass of gas.^° LIGs are 
generally regarded as recent galaxy mergers in which much of the gas of the colliding 
objects has fallen into a comm on ce nter (typically less than 1 kpc in extent), triggering 



a huge starburst phenomenon [344| . There is evidence for the existence of even more 



The Southern Site of PAO has been christened Pampa AmariUa. Recall that it has an aperture 
A Ri 7000 km^ sr for showers with incident zenith angle less than 60°. 
^° For stars to form, a large mass fraction at high density, dM/d{\ogn), is required. The Milky Way 
has an order of magnitude more gas at less than 300 cm~^ than what it has at 10** cm~^, contrary 
to LI Gs, which have nearly equal mass per decade of density between 10^ and 10* cm^"^ (see, e.g., 
Ref. [33, Im^) 
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extreme enviroments within LIGs (see, e.g., |345l |): These, larger than giant molecular 
clouds but with densities found only in small cloud cores, appear to be the most out- 
standing star formation places in the universe. They are well traced by HCN emission, 
i.e., they produce a substanti al fractio n of the whole HCN emission observed for the 
whole gXy H. see also H EHU- The CR enhancement factor in these small 
but massive regions can well exceed the average value for the galaxy. In A rp 220, for 
instance, two such regions were discovered to contain about 2 x 10^ Mq [345^ . If the CR 
enhancement in these regions is significantly larger than the starburst average, these 
extreme environments could be the main origin for any 7-ray emission observed from 
this galaxy j352!]. 



2. Propagation and further studies 



Using the PSCz catalogue [353|, Smialkowski et al. |354l | constructed all-sky maps 
of UHE proton intensities plausibly originating in LIGs, taking into account effects of 
particle propagation through the extragalactic medium and the possible influence of 
the regular galactic magnetic field. The PSCz catalogue consists of almost 15000 IR 
galaxies with known redshifts, covering 84% of the sky. Finding correlations with such 
an overdistributed sample might be thought as a risky business. There are, however, 
some phenomenological reasons — apart of LIGs being super-starbursts — to search for 
a possible UHECR origin in LIGs. 

Arp 299, one of the the brightest infrared source within 70 Mpc and a system of 
colliding galaxies showing intense starburst, appeared earlier as VV 118 in the list of 



candidates for the AG AS A triplet presented by in Ref. Il37 . Il40]. Indeed, for 
Arp 299, even when a hidden AGN was observed, it can not account for the whole 
FIR luminosity [355^ and a strong starburst activity is required t o exp lain it. It might 
be considered as a prime candidate for the origin of the triplet j354l | . if such is not a 
statistical fluctuation. 

Equatorial maps of the expected intensity of the UHE protons originating in LIGs for 
the energy region 40-80 EeV with sky coverage a nd d eclination dependent exposure for 
the AGASA experiment were presented in Ref. jsS^. Expected proton intensities for 
40-80 EeV appear to show good correlation with the distribution of the experimental 
events from AGASA, especially, from the region of the sky near Arp 299 and the 
AGASA triplet of events (RA^ 170°, 6 ^ 60°). No correlation of the data events 
above 80 EeV with the expected CR inten sitie s was reported. 

The latter result was confirmed in Ref. |203l | , where application of the Kolmogorov- 
Smirnov (KS) to the 11 highest AGASA events above 10^° eV yields a KS probability 
of < 0.5%, rejecting the possible association at > 99.5% significance level. This was 
used to argue that the existing CR events above 10^° eV do not owe their origin to 
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long burst GRBs, rapidly rotating magnetars, or any other events associated with core 
collapse supernovae (although it would yet be too early to completely rule out such a 
possibility, based only in a 2(T deviation). 

A more sensitive approach to study a possible LIG origin of UHECRs, perhaps, is 
not to look for possible correlations between the whole PSCz calatog and UHECRs, 
but rather to select a priori which LIGs are most likely to be detectable by their 
possible UHECR emission. There are several LIGs for which reasonable values of CR 
enhancements, comparable to, or lower than, the ratio between their SFR and the Milky 
Way's, can provide a 7-ray flux above GLAST sensitivity, and if the CR spectrum is 
sufficiently hard, also above the sensitivities of the new Cerenkov telescopes. These 
LIGs are then most likely to appear as new point-like 7-ray sources. Even when it 
is natural to expect that a LIG will emit 7-rays, only the more gaseous, nearby , and 



CR-enhanced galaxies are the ones which cou ld be detected as 7-ray point sources |852 

Out of the HCN just presented in Refs. |35d IsSlj ^^ and the larger Pico dos Dias 
Survey (PDS, |358l |)^^. the most likely 7-r av so urces are listed in Table HI together with 
EGRET fluxes upper limits, obtained in (352| . In our opinion, this group, as well as 
the HCN and the PDS should be separately taken into account when searching for 



•^^ Following |20,'j | . the core colla pse s upernova rate per galaxy, summed over all galaxy types, is 



5*311 ~ 0.011 SN per galaxy-year |205| . which yields a volume averaged rate of supernovae of S'sn'^g ~ 
2.2 X 10^'' SN/Mpc"^-year. The supe rnova rate for galaxies with far infrared emission is S^^ = 2.5 x 
IO^^'Lio SN/(FIR galaxy)-year |35(il | . Integrating this supernova rate over the FIR galaxy luminosity 
function yields 4*L'"firg ~ 0.7 x 10-'^Xobsc(i?o''''/300) SN/Mpc^-year. i^obsc, the correction for 
supernovae missed in the existing optical and near inf rared supernova detection surveys, might be 
as large as 10, and probably is at least as large as 3 |356j . Thus, the luminous infrared galaxies 
contribute at least 28% (i^obsc — 1) of the total supernova rate, with a total space density only 25% 
that of all normal galaxies. Furthermore, since the brightest infrared galaxies (Lfir > 10^^ Lq) 
dominate the cont ribution from all FIR galaxies, and these are quite rare, with space density ~ 4 x 
10~^ Mpc~^ correlations of UHECR arrival directions with the sky positions of the luminous 

IRAS galaxies offers a promising opportunity to test the hypothesis that UHECR acceleration has 
something to do with core collapse supernovae, as is implied by the GRB shock and magnetar 
unipolar inductor models for the acceleration sites. 

This survey is a systematic observation (essentially, all galaxies with strong CO and IR emission 
were chosen for HCN survey observations) of 53 IR-bright galaxies, including 20 LIGs with LpiR > 
IO-'^-'^Lq, 7 with LpiR > IO^^Lqj and more than a dozen of the nearest normal spiral galaxies. It 
also includes a literature compilation of data for another 9 IR-bright objects. This is the largest 
and most sensitive HCN survey (and thus of dense interstellar mass) of galaxies to date. 
The PDS survey consists of relatively nearby and luminous starbursts galaxies selected in the FIR. 
PDS galaxies have a lower mean IR luminosity log(LiR/L0) — 10.3 ± 0.5, redshifts smaller than 
0.1, and form a complete sample limited in flux in the FIR at 2 x 10~^° erg cm~^ s^^. 
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TABLE I: Powerful local LIGs (all with interferometric measurements) likely to be detected 
by GLAST, with EGRET upper limits. The luminosity distance in an standard universe 
{Dl = c/Hoqo^[l -qo + Qqz + (go - l)(2go2 + 1)^/^], ^0 (~75 km s"^ Mpc'^) is the Hubble 
parameter, qq (~0.5) is the deceleration parameter, and z is the redshift), the central-sphere 
radius from which the line emission was detected, the FIR luminosity and gas mass are also 
given. The minimum average value of CR enhancement k for which the 7-ray flux above 100 
MeV is above 2.4 xlO~^ photons cm~^ s~^, i.e. GLAST sensitivity is given. The smaller the 
value of (k) , the higher the possibility for these Galaxies to appear as 7-ray sources. 



Name Dl R \og{L^i^/\.Q) logM(H2/M0) {k) ^^Sm^cV 

[Mpc] [pc] [10~^ photons cm"^ s^-*^] 
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NGC1068 


15 


0.10 


10.74 


9.46 


78 


<3.6 


NGC2146 


20 


0.33 


10.78 


9.43 


149 


<9.7 


NGC4038/9 


22 


0.49 


10.65 


9.07 


412 


<3.7 


NGC520 


29 


0.38 


10.58 


9.47 


285 


<4.6 


IC694 


41 


0.30 


11.41 


9.59 


432 


<2.2 


Zw049.057 


52 


0.40 


10.95 


9.67 


578 


<6.9 


NGC1614 


64 


0.60 


11.25 


9.78 


680 


<5.0 


NGC7469 


65 


0.82 


11.26 


9.89 


544 


<3.2 


NGC828 


72 


0.92 


11.03 


10.09 


421 


<6.1 


Arp220 


72 


0.16 


11.91 


10.43 


193 


<6.1 


VV114 


80 


0.93 


11.35 


10.03 


597 


<3.9 


Arpl93 


94 


0.25 


11.34 


10.22 


532 


<5.2 


NGC6240 


98 


1.65 


11.52 


10.03 


896 


<6.4 


Mrk273 


152 


0.13 


11.85 


10.33 


1081 


<2.3 


IRAS 17208-0014 


173 


0.75 


12.13 


10.67 


640 


<7.5 


VIIZw31 


217 


1.27 


11.66 


10.70 


940 


<3.2 



possible UHECR correlations with new sets of data. 
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A stacking procedure with EGRET data is to be reported by Cillis et al. |35 
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H. Gamma-ray bursts 



1. Basic phenomenology 

Gamma ray bursts (GRBs) are flashes of liigli energy radiation tliat can be brigfiter, 
during tlieir brief existence, than any other gamma ray source in the sky. The bursts 
present an amazing variety of temporal profile s, sp ectra, and timescales that have puz- 
zled astrophysicists for almost three decades (86fll |. In recent years, our observational 
insight of this phenomenon has been dramatically improved by the huge amount of data 
collected by the Burst and Transient Source Experiment (BATSE): several thousand 
GRB observations were obtained. New breakthrough results are the expected outcome 
of HETE-2 and Swift. 

The temporal distribution of the bursts is one of the most striking signatures of the 
GRB phenomenon. There are at least four classes of distributions, from single-peaked 
bursts, including the fast rise and expon ential de caying (FREDs) and their inverse 



(anti-FREDs), to chaotic structures (e.g. |361l . l362| ). Burst timescales go through the 
30 ms scale to hundreds of seconds. 

The GRB photon spectrum is well fitted in the BATSE detectors range, 20 keV to 
2 MeV H, 

by a combination of two power-laws, dn^/de^ (x e^" {a is the flux 



density spectral index, Fj, oc with different values of a at low and high energy 
j363] . Here, dn^/de^ is the number of photons per unit photon energy. The break 
energy (where a changes) in the observer frame is typically e^b ~ 1 MeV, with a ^ 
at energies below the break and a ^ — 1 above the bre ak. In sev eral cases, the spectrum 
has been observed to extend to energies > 100 MeV 



The angular distribution of these bursts is isotropic, and the paucity of compara tively 



faint bursts implies that we are seeing to near the edge of the source population [365 
Both effects, isotropy and non-homogeneity in the distribution, strongly suggest a cos- 
mological origin of the phenomenon, confirmed by the detection of afterglows, delayed 
low energy emission of GRBs that allowed the measurement of the distance to the 
burst via a redshift determination of several GRB host-galaxies (e.g. 36d 367| ). 



2. The fireball model 



The most popular interpretation of the GRB-phenomenology is that the observable 
effects are due to the dissipation of the kinetic energy of a relativist i c expandi r ig plasma 
wind, 



a "fireball", whose primal cause is not yet known [368 . 
see [375^ for a detailed review). The rapid rise time and short duration 



1 ms 

of the burst imply that the sources are compact, with a linear scale comparable to a 
light-ms, ro ~ 10^ cm. If the sources are so distant, the energy necessary to produce 
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the observed events by an intrinsic mechanism is astonishing: about 10^^ erg of 7 
rays must be released in less than 1 second. Compactness and high 7-ray luminosity 
implied by cosmological distances result in a very high optical depth to pair creation, 
since the energy of observed 7-ray photons is above the threshold for pair production. 
The number density of photons at the source is 

= A-Krlcn-ye, (58) 

where e ~ IMeV is the characteristic photon energy. Using tq ~ 10'' cm, the optical 
depth for pair production at the source is 

r^^ ~ r^n^ax ~ ^^-^ ~ 10^^ . (59) 
47rroce 

The high optical depth creates the fireball: a thermal plasma of photons, electrons, 
and positrons. The radiation pressure on the optically thick source drives relativistic 
expansion, converting internal energy into the kinetic energy of the inflating shell |369l 
|370|. As the source expands, the optical depth is reduced. If the source expands with 
a Lorentz factor F, the energy of photons in the source frame is smaller by a factor F 
compared to that in the observer frame, and most photons may therefore be below the 
pair production threshold. 

Baryonic pollution in this expanding flow can trap the radiation until m ost of th e 
initial energy has gone into bulk motion with Lorentz factors of F > 10^ — 10^ 375 . 376l . 



The kinetic energy, however, can be partially converted into heat when the shell collides 
with the interstellar medium or when shocks within the expanding source collide with 
one another. The randomized energy can then be radiated by synchrotron radiation 
and inverse Compton scattering yielding non-thermal bursts with timescales of seconds, 
at large radius r = > lO^^cm, beyond the Thompson sphere. Relativistic random 
motions are likely to give rise to a turbulent build up of magnetic fields, and therefore 
to Fermi acceleration of charged particles. 

Coburn and Boggs 1377^ recently reported the detection of polarization, a particu- 
lar orientation of the electric- field vector, in the 7-rays observed from a burst. The 
radiation released through synchrotron emission is highly polarized, unlike in other 
previously suggested mechanisms such as thermal emission or energy loss by relativis- 
tic electrons in intense radiation fields. Thus, polarization in the 7-rays from a burst 
provides direct evide nce i n support of synchrotron emission as the mechanism of 7-ray 
production (see also |378l |). 

3. Fermi acceleration in dissipative wind models of GRBs 

Following Hillas' criterion, the Larmor radius should be smaller than the largest 
scale Zgrb over which the magnetic field fluctuates, since otherwise Fermi acceleration 
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will not be efficient. One may estimate /grb as follows. The comoving time, i.e., 
the time measured in the fireball rest frame, is t = r/Tc. Hence, the plasma wind 
properties fluctuate over comoving scale length up to Iqkb ~ ''^/r, because regions 
separated by a comoving distance larger than r/T are causally disconnected. Moreover, 
the internal energy is decreasing because of the expansion and thus it is available for 
proton acceleration (as well as for 7-ra y pro duction) only over a comoving time t. The 



typical acceleration time scale is then 248 1 



where j3c is the Alfven velocity. In the GRB scenario /5 ~ 1, so Eq. (jUn|) sets a 
lower limit on the required comoving magnetic field strength, and the Larmor radius 
r-^ = E'/eB = E/TeB, where E' = E/T is the proton energy measured in the fireball 
frame. 

Th is condition sets a lower limit for the required comoving magnetic field strength 



> 0.151^00^20^51 , (61) 



where E = 10^''£'2oeV, F = SOOFsoo, L = lO^^Lsierg s~^ is the wind luminosity, and 
Bg p is the equipartition field, i.e. a field with comoving energy density similar to that 
associated with the random energy of the baryons. 

The dominant energy loss process in this case is synchrotron cooling. Therefore, the 
condition that the synchrotron loss time of Eq. pip be smaller than the acceleration 
time sets the upper limit on the magnetic field strength 

5 < 3 X 10'Tl,,E^,^G. (62) 

Since the equipartition field is inversely proportional to the radius r, this condition 
may be satisfied simultaneously with (|?)T|l provided that the dissipation radius is large 
enough, i.e. 

r, > 10''T^,lEl,cm. (63) 

The high energy protons lose energy also in interaction with the wind photons (mainly 
through pion production). It can be show n, h owever, that this energy loss is less 



important than the synchrotron energy loss |248 . 

A dissipative ultra-relativistic wind, with luminosity and variability time implied by 
GRB observations, satisfies the constraints necessary to accelerate protons to energy 
> 10^° eV, provided that T > 100, and the magnetic field is close to equipartition 
with electrons. We stress that the latter must be satisfied to account for both 7-ray 



emission and afterglow observations |376| . At this stage, it is worthwhile to point out 
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that for the acceleration process at shocks with large T the particle distributions are 
extremely anisotropic in shock, with the particle angular distribution opening angles 
~ in the upstream plasma rest frame. Therefore, when transmitted downstream 
the sh ock p articles have a limitted chance to be scattered efficiently to re-eneter the 
shock [2221 • However, in this particular case, the energy gain by any "successful" CR 
can be comparable to its original energy, i.e., {AE)/E ~ 1. 



4. UHECRs and GRBs: connections 



In the GRB model for UHECR production described above^^, the high energy CRs 
are protons accel erated by Fermi's mechanism in sources that are distributed through- 
out the universe 248l l387| . It is therefore possible to compare the UHECR spectrum 
with the prediction from a homogeneous cosmological distribution of sources, each gen- 
erating a power law differential spectrum oc of high energy protons as typically 
expected from Fermi acceleration. Under the assumption that the GRB rate evolution 
is sim ilar to the star-formation rate evolution, the local GRB rate is ~ 0.5 Gpc~^ yr~^ 
|388j |. implying a local 7-ray energy generation rate of ~ 10^"^ erg Mpc~^ yj,-i_36 

The energy observed in 7-rays reflects the fireball energy in accelerated electrons. If 
accelerated electrons and protons (as indicated by afterglow observations \38^ ) carry 
similar energy, then the GRB production rate of high energy protons is 

el{dhp/dep),=o ~ W erg Mpc-=^ (64) 

The generation rate (Eq. of high energy protons is remarkably similar to that re- 
quired to account for the fiux of > 10^^ eV CRs, whereas in this model, the suppre ssion 
of model fiux above 10^^ '' eV is due to the GZK cutoff. Stecker and Scully |j39^ ISQll ] 
have raised doubts on the possibility of this generating a very strong cutoff at the 
highest CR energies, since if the GRB redshift distribution follows that of the star 
formation rate in the universe, a rate which is higher at larger redshift, most of the 
GRBs would be just too far and CR with energies above 3 x 10^^ eV would be strongly 
attenuated by the CMB. For a HiRes-shape spectrum, a common origin between GRBs 



An alternative model for the GRB phenomenon has been recently put forward [38(fi- In such a 
model, the GRB explosion occurs when a massive star collapses into an electromagnetic blac k hol e 
which readily discharges due to the annihilation of pairs produced by vacuum polarization [sSlj . 
The model e xpla i ns the tim e variability, the spectra, and the GRB afterglows to a very high level 
of accuracy I383I l384| . Interestingly, within this set up one is ab le to acce lerat e the ionized 



hydrogen atoms sorrounding the death star to ultra high energies [sssf (see also 386] ) . 
The local (z = 0) energy production rate in 7-rays by GRBs is roughly given by the product of the 
characteristic GRB 7-ray energy, E fa 2 x 10^^ erg, and the local GRB rate. 
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and ultrahigh energy CRs j392j | is favored. For ap praisals of this and other general 
criticisms made to the GRB-UHECR connection see 393 . 39^. 

Two of the highest energy CRs come from directions that are within the error boxes 
of two rema rkable GRBs detected by BATSE with a delay of (9(10) months after the 
bursts |395l | . However, a rigorous analysis shows no correlation between the arrival 
direction of UHECRs and GRBs from the third BATSE catalog |,396.] . No correlations 
were found either between a pre-CGRO burst catalog and the Haverah Park shower set 
that covered approximately the same period of time. These analysis, however, could 
have been distorted by the angular resolution {A6 ~ 3°) of the GRB measurements. A 
sensitive anisotropy analysis between UHECRs and GRBs will be possible in the near 
future, using PAO, HETE-2 and Swift. Preliminary results (if one assumes that GRBs 
are most likely to happen in infrared luminous galaxies) do not seem to indicate any 
strong correlation (see above Section Hi Gj) . 



5. A GRB origin for CRs below the ankle? 

Wick, Dermer and Atoyan |397| have recently proposed a model for the origin of 
all CRs above ~ 10^^ eV/nucleon. In this model, GRBs are assumed to inject CR 
protons and ions into the interstellar medium of star-forming galaxies -including ours- 
with a power-law spectrum extending to a maximum energy ~ 10^° eV. In addition 
to the more energetic, extragalactic spectrum of CR, the CR spectrum near the knee 
was also shown to be plausibly fitted with CRs trapped in the Galactic halo that 
were accelerated and injected by an earlier Galactic GRB.^^ For power-law CR proton 
injection spectra oc E""^'"^ and low and high-energy cutoffs, normalization to the local 
time- and space-averaged GRB luminosity density implies that if this model is correct, 
the nonthermal content in GRB blast waves is hadronically dominated by a factor ^ 60- 
200, limited in its upper value by energetic and spectral considerations. Neutrinos to 
be detected in kilometer-scal e neu trino detectors such as IceCube (See Sec.lIV|) provide 
a clean signal of this model |400| . 

The last GRB in the Galaxy has been also proposed as the possible progenitor of 



In addition, dispersion of magnetic fields in the intergalactic medium can make the number of 
UHECR-contributing GRBs to grow above the burst rate within the GZK sphere. The latter, 
within 100 Mpc from Earth, is in the range of 10^^ to lO^'^ yr~^. Assuming a dispersion tiemscale, 
~ lOjyr, the number of sources contributing to the flux at any given time may be as large as 

Arguably, from the evidence for beaming and the association of GRBs with star-forming galaxies 
like the Milky W ay, GRB events are estimated to occur once every 10'' — 10^ yrs in the inner 
Galaxy IsoilsQfll . 
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the CR anisotropy observed in the direction of the GC |40]| . In order to estimate the 
remaining traces of any CR activity produced by the last GRB in the galaxy, one has 
to take into account several considerations: 
(t) The UHECRs escaping the GRB fireball 



No{E > eV) ~ 10"' No{E > E^-.^) (65) 

are mostly neutrons, because protons a re ca ptive in the magnetic field and suffer ex- 
tensive adiabatic losses on the way out 402] -^^ Some of these neutrons will decay into 



protons within the GC thin disk-like (r ~ 3 kpc) region of high interestellar medium 
density and high star formation rate. The population of secondary protons would then 
be captured by the strong 5-field near the GC, attaining diffusion with a residence 
time scale of about T ~ 10^ yr. At the end of this time, about 1/300 protons are able 
to avoid leakeage. The trapped protons, 

Nt{E > 10^^ eV) = No{E > 10^^ eV) (^1 - e""^) /300 , (66) 

can then be turned back into neutrons by interaction with nuclei in the interstellar 
medium with probability of 5 x 10"^. Here, m„ and r are the neutron mass and lifetime, 
respectively. 

(ii) The formation of the n ^ p reservoir depends on the GRB rate in the inner Galaxy 
(taken here as 10^ yr~^) times the probability that a GRB jet points more o less along 
the direction of the GP. The latter is estimated to be about 50%. 
(in) The total CR production by a single GRB is ~ 10^^ erg |403l |. 

Now, recalling that the CR excess observed by AGASA represents a luminosity 
of particles beyond 10^^ eV of about 4 x 10^° erg/s [59], straightforward calculation 
shows that the observed anisotropy in the direction of the GC can be easily fitted 
by the neu trons produced in the GRB reservoir, that ultimately travel unscathed to 



Earth 401 



III. FULL THROTTLE: UHECR GENERATION BEYOND THE STAN- 
DARD LORE 

The astrophysical models discussed present difficulties in providing a completely 
satisfactory explanation of the super-GZK events, if there are any. This may simply 
reflect our present lack of statistics, our present ignorance of the true conditions of 
processes in some highly energetic regions of the universe, or, perhaps, may imply 



We remind the reader that the differential injection spectrum of GRBs cx E 
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that exotic mechanisms are at play. Physics from the most favored theories beyond 
the standard model (SM) hke string/M theory, supersymmetry (SUSY), grand unified 
theories (GUTs), and TeV-scale gravity have been invoked to explain the possible flux 
above the GZK energy limit. This review is not mainly concerned with beyond-SM 
scenarios (for more comprehensive surveys see e.g. [H 0, 4Q4, 405 . 1406.] ). but for the 
sake of completeness, we provide here a brief account of some of the most relevant 
exotic explanations. 

The most economical among hybrid proposals involves a familiar extension of the SM, 
namely, neutrino masses. It was noted many years ago that z/'s arriving at Earth from 
cosmologically distant sources hav e an annihilation probability on t he r elic neutrino 
background of roughly 3 hQ^% |407| . Inspired on this analysis, Weiler |408| and Fargion 



et al. |409| noted that neutrinos within a few Z widths of the right energy, 

ERz = ^ = ^ f — 1 X 10^1 eV , (67) 
to annihilate with the relic neutrinos at the Z-pole with large cross section. 



(flann)^ = I ^ <y^Us) = 2 TT V2Gf ~ 40.4 ub , (6^ 

z 



may produce a "local" flux of nucleons and photons.^'' Remarkably, the energy of the 
neutrino annihilating at the peak of the Z-pole has to be well above the GZK limit. The 
mean energies of the ~ 2 nucleons and ~ 20 7-rays in each process can be estimated by 
distributing the resonant energy among the mean multiplicity of 30 secondaries. The 
proton energy is given by 

^ J^L ^ 1.3 10^° eV, (69) 



whereas the 7-ray energy is given by 



Ml / eV ^ ,n 



"3 



The latter is a factor of 2 smaller to account for the photon origin in two body tt" 
decay. This implies that the highly boosted decay products of the Z could be observed 
as super-GZK primaries. However, to reproduce the observed spectrum, the Z- 
burst mechanism requires very luminous sources of extremely high energy neutrinos 



"0 Gf = 1.16639(1) X 10"^ GeV-2 is the Fermi coupling constant. 

Similarly, gravi-burst fragmentation jets can contribute to the super-GZK spectrum |41( 
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411 



412 



413|. The p resen t limits on these sources are near 
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throughout the universe 

the threshold of sensitivity to the required flux 

In 1931, Georges Lemaitre 415l | - a forerunner of the Big Bang hypothesis - intro- 
duced the idea that the entire material fllling the universe, as well as the universe's 
expansion, originated in the super-radioactive disintegration of a "Primeval Atom", 
which progressively decayed into atoms of smaller and smaller atomic weight. The CRs 
were introduced in this picture as the energetic particles emitted in intermediate stages 
of the decay-chain. Echoing Lemaitre, in the so-called "top-down models" , charged and 
neutral prin aaries arise in the quantum mechanica l decay of sup e rma s sive elenaentar y 

X particles 113, H £3, £3, Q H, H H, H H, H, H H H H- 

To maintain an appreciable decay rate today, it is necessary to tune the X lifetime to 
be longer (but not too much longer) than the age of the universe, or else "store" short- 
lived X particles in topological vestiges of early universe phase transitions (such as 
magnetic monopoles, superco nducting cosm ic strings, vortons, cosmic ne cklaces, etc.). 



Discrete gauged symmetries |43lL I432L l433| or hidden sectors |434 l435| are generally 
introduced to stabilize the X particles. Higher dimensional operators, wormholes, 
and instantons ar e then inv oked to break the new symmetry super-softly to maintain 
the long lifetime 423 . l424l | (coUissional annihilation has been considered too |43fil | ) . 
Arguably, these metastable super-heavy relics (MSRs) may constitute (a fraction of) 
the dark matter in galactic haloes. 

Of course, the precise decay modes of the X's and the detailed dynamics of the 
flrst generation of secondaries depend on the exact nature of the X particles under 
consideration. However, in minimal extensions of the SM, where there are no new 
mass scales between Msusy ~ 1 TeV and m^, the squark and sleptons would be- 
have like their corresponding supersymmetric partners, enabling one to infer from the 
"known" evolution of quarks and leptons the gross features of the X particle decay: 
the strongly interacting quarks would fragment in to ie t s of hadro ns containing mainly 



437, 



438L |439|.^^ This implies that 



pions together with a 3% admixture of nucleons 
the injection spectrum is a rather hard fragmentation-type shape (with an upper limit 
usually fixed by the GUT scale) and dominated by 7-rays and neutrinos produced via 
pion decay. Therefore, the pho ton/proton ratio can be used as a diagnostic tool in 
determining the CR origin (442j |. In light of the mounting evidence that UHECRs are 
not 7-rays, one may try to force a proton dominance at ultrahigh energies by postulat- 
ing efficient absorption of the dominant ultrahigh energy photon flux on the universal 
and/or galactic radio background. However, the neutrino flux accompanying a nor- 
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In light of this, the s ensitivity of future CR experiments to test the SUSY parameter space has been 
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Mil- 



estimated 

^'^ Of course, this is not the case in traditional scenarios (for an exception see |443| 1 of decaying massive 
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malized proton flux is ine vitab ly increased to a level where it should be within reach 
of operating experiments 
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It is clear that because of the wide variety of top down models the ratio of the 
volume density of the X-particle to its decay time is model dependent. However, if a 
top down scenario is to explain the origin of UHECRs, the injection spectrum should 
be normalized to account for the super-GZK events witho ut vi olating any observational 
flux measurements or limits at higher or lower energies 445[. In particular, neutrino 



and 7-ray fluxes depend on the energy released integrated over redshift, and thus on 
the specific top down model. Note that the electromagnetic energy injected into the 
Universe above the pair production threshold on the CMB is recycled into a generic 
cascade spectrum below this threshold on a short time scale compared with the Hubble 
time. Therefore, it can have several potential observable effects, such as modified light 
element abundances due to ^He photodisintegratipn, or induce spectral distortions of 
universal 7-ray and neutrino backgrounds 446l l447l | . Additionally, measurements of the 



diffuse GeV 7-ray flux |448l |. to which the generic cascade spectrum would contribute 



directly, limit sign i ficantly t he parameter space in which X's can generate the fiux of 
the UHECRs 



449, 



especially if there is already a significant contribution to 



this background from conventional sources such as unresolved 7-ray blazars. Recently, 
a possible lower extragalactic c ontributio n to the diffuse 7-ray background measured by 
EGRET has been pointed out llii,^. The ~ 50% smaller EGRET flux practically 
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rules out extragalactic top down and Z-burst scenarios 

If MSRs are the progenitors of the observable UHECRs, then the flux will be domi- 
nated by the decay or annihilation products of X's in the Galactic halo, i.e., by sources 
at distances smaller than all relevant interaction lengths. A clean signal of this sce- 
nario is the pred icted anisotropy due to the non-central position of the Sun in our 
galaxy |455|, |456|.^'' Although it was noted that the predict ed a nisotropy is consistent 
with AGASA and Haverah Park data recent analyses |458l . L459] of SUGAR data 
excludes the MSR hypothesis at the 5a level if all events above 10^^'^ eV are due to 
metastable X clustered in the halo. (For the extreme case where the population of 
MSRs is responsible only for CRs with energies > 10^^'^ eV, the annihilation scenario 
is disfavored at least at the 99% CL, whereas decaying MSRs still have a probability 
~ 10% to reproduce SUGAR data.) 

A more exotic explanation postulates that the X's themselves constitute the pri- 
maries: magnetic monopoles easily pick up energy from the magnetic flelds permeating 



dark matter in the Galactic halo, which due to the lack of absorption, predict compositions directly 
given by the fragmentation function, i.e., domination by 7-rays. 

Additionally, the arrival directions of UHECRs can show significant deviations from a random 
distribution due to an anisotropic spectrum of relics in the early universe j453|. 
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the universe and can traverse unscathed through the p rimev al radiation, providing an 
interesting candidate to generate extensive air showers 46Cll |. In particular, a baryonic 
monopole encountering the atmosphere will diffuse like a proton, producing a compos- 



ite heavy-particle-like cascade after the 
muons among all the charged particles 
a poorly understood super-GZK event 



irst interaction j46ll | with a great number of 
46^ . Although this feature was observed in 
46l |464], it 

explanation for the UHECR data sample would be in terms of magnetic monopoles 
alone. Moreover, any confirmed directional pairing of events would appear difficult to 
achieve with the monopole hypothesis. 

A novel beyond-SM-model proposal to break the GZK barrier is to assume that 
UHECRs are no t known particles but rather a new species, generally referred to as 
the uhecron, U 



465 



466l . 14671] . The meager information we have about super-GZK 



particles allows a naive description of the properties of the U. The muonic content in 
the atmospheric cascades suggests f/'s should interact strongly. At the same time, if 
U's are produced at cosmological distances, they must be stable, or at least remarkably 
long lived, with mean-lifetime r > 10^ {mu/3 GeV) (rf/Gpc) s, where d is the distance 
to the source and mjj, the uhecron's mass. Additionally, since the threshold energy 
increases linearly with mu, to avoid photopion production on the CMB mr r r*, 1.5 GeV. 



In recent years, direct searches of supersymmetric hadrons j468L 146 9L l47flL l47ll | have 
severely eroded the attractivene ss of the U scenario. However, adequate fine-tunings 
leave a small window still open 

On a similar track, it was recently put forward 
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that strangelets (stable lumps 
of quark matter with roughly equal numbers of up, down, and strange quarks) can 
circumvent the acceleration problem in a natural way (due to a high mass and charged, 
but low charge-to-mass ratio) and move the expected cutoff to much higher energies. 

Another possibility in which super-GZK CRs can reach us from very distant 
sourc es ma v arise out of photons that mix with light axions in extragalactic magnetic 
fields 1475^ . These axions would be sufficiently weakly coupled to travel large distances 
unhindered through space, and so they can convert back into high energy photons 

ocal 
In 



close to the Earth. An even more radical proposal postulates a tiny violation of 
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Lorentz invariance, such that some processes become kinematically forbidden 
particular, photon-photon pair production and photopion production may be affected 
by Lorentz invariance violation. Hence, the absence of the GZK-cutoff would result 
from the fact that the threshold for photopion production "disappears" and the pro- 
cess becomes kinematically not allowed. This implies that future PAO observations of 



See also |476j | for a SUSY inspired related scenario. 
Existing limit s on the violation of Lorentz invariance come from multi-TeV 7-ray observations of 
Mrk 501 l47i. 
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faraway sources could provide constraints on, or even a measurement of, the violation 
of Lorentz symmetry, yiel ding e ssen t ial insights into the nature of grayity- in duced wave 
dispersion in the vacuum jili H, H, mi H, Esi, isl H, H . 

In summary, future UHECR data may not only provide clues to the particle sources, 
but could enhance our understanding of fundamental particle physics. We are enterin g 
this new High Energy Physics era with the Pierre Auger Observatory [isi H, liill 



IV. MINORITY REPORT: NEUTRINO SHOWERS 

Extraterrestrial neutrinos provide a unique window to probe the deepest reaches 
of stars, quasars, and exotic structures in the cosmos. In contrast to all other SM 
particles, they can escape from dense astrophysical environments and propagate to the 
Earth unscathed, arriving aligned with their source. Even at the highest energies the 
uu annihilation mean free path on the cosmic neutrino background, 

A. = (n, a,p)-^ ^ 4 X lO^^ cm , (71) 
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is somewhat above the present size of the horizon (recall that /Jq ^ ~ 10^^ cm) 

In addition, the fluxes of ultra-high energy cosmic neutrinos offer clues to the prop- 
erties of neutrinos themselves and provide an important probe of new ideas in particle 
physics. Namely, their known interactions are so weak that new physics may easily 
alter neutrino properties. This is especially relevant for extraterrestrial neutrinos with 
energies ^10^ GeV that interact with nucleons in the Earth atmosphere with center- 
of-mass energies above the electroweak scale, where the SM is expected to be modified 
by new physics. 

In summary, ultra high energy extraterrestrial neutrinos are inextricably linked with 
the physics processes discussed throughout this review. This chapter will therefore 
focus on the main features of high energy neutrino astronomv. For comprehensive 



reviews on cosmic neutrinos the reader is referred to L49 



494, 
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A. Bounds on the energy spectrum 

The basic operational system of neutrino telescopes is an array of strings with photo- 
multiplier tubes (PMTs) distributed throughout a natural Cerenkov medium such as 
ice or water. The largest pilot experiments (~ 0.1 km in size) are: the now defunct 
DUMAND (Deep Underwater Muon and Neutrino Detector) ex perim ent |497l | in the 
deep sea near Hawaii, the underwater experiment in Lak e Baikal |49q , and AMANDA 
(Antarctic Muon And Neutrino Detector Array ) \A9^ in the South Pole ice. Next 
generation neutrino telescopes aim towards an active volume in the range of 1 km^ 
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of water. Projects under construction or in the proposal stage are: two deep sea 
experiments in the Mediterranean, the Fre nch A NTARES (Astronomy with a Neutrino 
telescope Abyss enviro nme nt RESearch) SOOl l and NESTOR (Neutrino Experiment 
SouthwesT Of GReece |50l|), and IceCube [SoJ, a scaled up version of the AMANDA 



detector. 

The traditional technique to observe cosmic neutrinos is to look for muons (along 
with a visible hadronic shower if the v is of sufficient energy) generated via charged- 
current interactions, {vy^^v^N — » (/i~,yU+)+ anything, in the rock below the detec- 
tor. The Cerenkov light emitted by these muons is picked up by the PMTs, and is 
used for track reconstruction. The muon energy threshold is typically in the range of 
10 — 100 GeV. However, besides the desired extraterrestrial neutrinos, below ~ 10^ GeV 
there is a s ignificant background of leptons (produced by CRs interacting in the at- 



mosphere) j503l . l504| and so the task of developing diagnostics for neutrino sources 
becomes com plica ted. With rising energy (> 10^ GeV) the three neutrino flavors can 
be identified j505|. 

The spectacular neutrino fireworks (in the 10 MeV r ange ) from supernova 1 987A 
constit ute t he only ext raga lactic source so f ar ob served 506 . l507l |. The Frejus 508 1. 
Baikal MACRO |m3, and AMANDA collaborations reported no excess of 

neutrinos above the expected atmospheric background, enabling significant limts to be 
set on the diffuse neutrino flux.^^ These limits are shown in Fig. El 

Upward going neutrinos with energies > 10^ GeV are typically blocked by the Earth. 
This shadowing severely restricts the high energy event rates in underground detectors. 
However, neutrinos may also induce extensive air showers, so current and future air 
shower experiments might also function as neutrino detectors. The neutrino interac- 
tion length is far larger than the Earth's atmospheric depth, which is only 0.36 km 
water equivalent (kmwe) even when traversed horizontally. Neutrinos therefore shower 
uniformly at all atmospheric depths. As a result, the most promising signal of neutrino- 
induced cascades are quasi-horizontal showers initiated deep in the atmosphere. For 
showers with large enough zenith angles, the likelihood of interaction is maximized and 
the background from hadronic cosmic rays is eliminated, since the latter shower high in 
the atmosphere. These z/-showers will appear as hadronic vertical showers, with large 
electromagnetic components, curved fronts (a radius of curvature of a few km), and 
signals well spread over time (on the order of microseconds). 



The Frejus experiment, located in an underground laboratory, measured the energy of muons pro- 
duced by neutrino interactions in the rock above the detector. The detector itself comprised a 
calorimeter made from a vertical sandwich of over 900 iron sl abs i nterspersed with flash chambers. 
Geiger tubes embedded in the structure provided the trigger |512l | . T he ex tensive air shower array 
on top of the Gran Sasso Laboratory (EAS-TOP) [sisl l and MACRO '514j used similar techniques. 
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524i | and ANITA (45 days running) •524'] corresponding to 1 event per decade. The 



runnmg) 

region between the falling dashed-dotted lines indicates the flux of atmospheric neutrinos. 



The event rate for quasi-horizontal deep showers from ultra-high energy neutrinos 



IS 
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N 



i,X 



^ / dEiNA-^ aiN^x{Ei)S{Ei) 



(72) 



where the sum is over all neutrino species i = z/g, z/g, i/^, z/^, z/^, u^, and all final states X. 
Na = 6.022 X 10^^ is Avogadro's number, and d^i/dEi is the source flux of neutrino 
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species i, a as usual denotes the cross section, and 8 is the exposure measured in 
cm^ w.e. sr time. The Fly's Eye and the AGASA Collaborations have searched for 
quasi-horizontal showers that are deeply -penetrating, with depth at shower maximum 
-^max > 2500 g/cm^ (see e.g. 525l l526| |). There is only 1 event that unambiguously 
passes this cut with 1.72 events expected from hadronic background, implying an upper 
bound of 3.5 events at 95%CL from neutrino fluxes. Note that if the number of events 
integrated over energy is bounded by 3.5, then the same limit is certainly applicable 
bin by bin in energy. Thus, using Eq. (ff^ one obtains 

/ dE, Na -j^ (T.^^xiE,) £iE,) < 3.5 , (73) 
J A dE, 

at 95% CL for some energy interval A. Here, the sum over X takes into account charge 
and neutral current processes. In a logarithmic interval A where a single power law 
approximation 

^^a,N^x{E.)S{E,)^Er (74) 
is valid, a straightforward calculation shows that 

<E)e^/^ dE,_ d<!>,, „ , ^ „ _ sinh5 



rkE)e^i^ ^E sinh5 , 

/ —^Ei—^ai^^x8 = {o,^^x8Eid^ildE,)^—^, (75) 



where 5 = (a + l)A/2 and {A) denotes the quantity A evaluated at the center of the 
logarithmic interval. The parameter a = 0.363 + /5 — 7, where 0.363 is the power law 
index of the SM neutrino cross sections ika and /3 and —7 are the power law indices 
(in the interval A) of the exposure and flux d^i/dEi, respectively. Since sinh 6/6 > 1, 
a conservative bound may be obtained from Eqs. (fTSj) and (fTSj) : 

Na Y,{a,N^x{Ei)) {S{E,)) {Ed^i/dE,} < 3.5/A . (76) 

i,X 

By taking A = 1 as a likely interval in which the single power law behavior is valid 
(this corresponds to one e-folding of energy) , and setting (Eid^i/dEi) = ^{E^d^^/dE^) 
($;^ = total neutrino flux) from Eq. (f7H|) it is straightforward to obtain 95%CL upper 
limits on the neutrino flux. The limits are shown in Fig. 121 

These bounds will be improved in the near future. In part icular, each site of PAO 
will observe ~ 15 km^we sr of target mass around 10^^ eV |52q . Moreover, the study of 
radio pulses from electromagnetic showers created by neutrino interactions in ice would 
provide an increase in the effective area up to lO^^km^. A prototype of this technique 
is the Radio Ice Cerenkov Experiment (RICE) [^29] . Similar concepts are used by 
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the Goldstone Lunar ultrahigh energy neutrino Experiment (GLUE),^^ the ANtartic 
Impulse Transient Array (ANITA), and the Fast On-orbit Recording of Transient 
Events (FORTE). 5° Existing limits on i/-fluxes from RICE, GLUE and FORTE, as 
well as the projected sensitivities of forthcoming experiments are collected in Fig. El 



B. Astronomy on Ice 

IceCube is, perhaps, the most promising route for neutrino detection |522l |. This 
telescope will consist of 80 kilometer-length strings, each instrumented with 60 10- 
inch photomultipliers spaced by 1.7 m. The deepest module is 2.4 km below the ice 
surface. The strings are arranged at the apexes of equilateral triangles 125 m on a 
side. The instrumented detector volume is a cubic kilometer. A surface air shower 
detector, IceTop, consisting of 160 Auger-style Cerenkov detectors deployed over 1 km^ 
above IceCube, augments the deep-ice component by providing a tool for calibration, 
background rejection and air-shower physics. Muons can be observed from 10^^ eV to 
10^^ eV. Cascades, generated by i^e, z^e, J^r, and Vr can be observed above 10^^ eV and 
reconstructed at energies somewhat above 10^^ eV. The angular resolution is ~ 0.7° at 
TeV energies. 

A vari etv o f neutrino-emitting-sources have been proposed, such as supernovae |53ll| 
(see also |532| for a detailed study o n dir ectional ity) , GRB fireballs 4od 533 ^534. 



microquasars j^yj. X-ray binaries |l69| . AGN |537j | . FRI radio galaxies etc. As 

an example of the IceCube potential, in what follows we briefly discuss its sensitivity 
to probe the neutron hypothesis of UHECRs via observation of the antineutrino beam 
n ^ p + e~ + z7e, expected from the Cygnus direction j^. To this end, we first estimate 
the background signal. As discussed in Sec. Ill At the TeV 7-ray flux. 



IT^ / F \ -l-9(±0.3stat±0.3sys) 

^ = 4.7(±2.l3tat ± l.Ssys) X 10-1=^ — =— phcm-'s-^TeV-^ , 

dE^ V =.ys; \^lTeVy ^ 

reported by HEGRA Collaboration |152l | in the vicinity of Cygnus 0B2 is likely due 
to hadronic processes. Since 7r°'s, vr+'s, and 7r"'s are produced in equal numbers. 



This experiment observes microwave C erenk ov pulses from electromagnetic showers induced by 
neutrinos interacting in the Moon's rim |52c| 

A balloon-borne payload that will circle the continent of Antartica at 35,000 meters, scanning 
the vast expanses of ice for telltale pulses of radio emission generated by the neutrino collisions 
[http : //www.ps .uci . edu^aLnita/) . 

The FORTE satellite records bursts of electromagnetic waves arising from ne ar th e Earth's surface 
in the radio frequency range 30 to 300 MHz with a dual, polarization antena |52ll |. 
Extension of these aperture is in the proposal stage 530l | . 
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we expect two photons, two i^e's, and four z/^'s per vr^. On average, the photons 
carry one-half of the energy of the pion, and the neutrinos carry one-quarter. For 
dFT^/dE^^ oc the energy-bins dE scale with these fractions, and we arrive at 

dF dF 



dE^' ' ' dE,, 

^(S„^.^^/4).4^(EJ, (78) 

for the fluxes at the source, wh ere tt denotes any one of the three pion charge-states. 
Terrestrial experiments (see e.g. jsS^) have shown that z/^ and Ut are maximally mixed 
with a mass-squared difference ~ 10~^eV^. This together with the known smallness 
of |(z/e|i/3)p, implies that the z/^'s will partition themselves equally between z/^'s and 
z/^'s on lengths large compared to the oscillation length Aqsc ~ 1-5 x 10^^ (E^/PeV) pc. 
Here z/3 ~ (z/^ + z/^) / is the third neutrino eigenstate. From these remarks, one finds 
a n early identical flux for each of the three neutrino flavors (j = e, fi, r), which is equal 

to liei] 

^(E,^E,/2)^2^(E,). (79) 

Although TeV neutrinos are copiously produced, be caus e they are weakly interacting 
the detection probability on Earth is tiny, about 10~^ |493^. In particular, the expected 
p event rate at IceCube associated with the unidentifled HEGRA source is < 1 yr~^ 
(D. Hooper, private communication). Such an event rate is even smaller than the 
atmospheric neutrino background. Moreover, existing limits on TeV 7-ray fluxes in 
this region of the sky are near the HEGRA sensitivity.^^ In light of this, we take as 
background the atmospheric neutrino event rate, and so Poisson statistics implies that 
a signal > 3.7 events is signiflcant at the 95% CL. 

Antineutrinos take only a very small part of the energy of the parent neutron, typi- 
cally ~ 10"'^. Hence, to estimate the event rate of TeV antineutrinos at IceCube, the 



For a year of running at IceCube the expected background from atmospheric neutrinos (with energy 
> 1 TeV) within 1° circle centered in the Cygnus direction (about 40° below the horizon) is < 1.5 
events. 

Specifically, the CAS A- MI A experiment observed the Cygnus region over a 5 year period (1990- 
1995) a nd th e collected data place a bound F^{E^ > 115 TeV) < 6.3 x 10"^^ phcm"^ s"\ at the 
90%CL |54ol |. Extrapolation down to lower energies, assuming a spectrum cx i?^^, leads to F^( E^ > 
ITeV) < 7.2 X 10^13 phcm-2g-i_ ^j^jg within a factor of 2 of the HEGRA measurements 
ImiI, F-y{E^ > ITeV) w 4.5 x 10"" g-i_ 
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relevant nucleus population at the source has an energy per nucleon En^psy ~ 1 PeV. 
Nuclei with Lorentz factor ~ 10^ are synthesized in a ll su pernovae. Hadronic inter- 
actions with the HII population (density < 30 cm~^ jlSS^ ) and photodisintegration 
processes provide the flux of PeV neutrons. In this energy regime, the target photons 
at photodisintegration threshold energies are in the ultraviolet, ~ 5 eV. This includes 
the entire emission spectrum of the O stars and about 60% of photons from B stars 
(with average temperature 28,000 K). From the photon emission rate Fjjy the number 
density riuv at the surface of a sphere of radius R from the core center is given by 

For the 0-star population, the photon emission rate in the Lyman region is found to 
be Fl ~ 10^^ photons s~^ [206J. The Lyman emission corresponds to 60% of the entire 
O star spectrum. Furthermore, as mentioned above, 60% of the B star spectrum is also 
active for photodisintegration in this energy regio n, and the B star population is about 



20 times greater than that of the O stars [206|. Now, from the H-R diagram |542 
one can infer that the energy luminosity of a B-star is about 0.1 that of an O star. 
Additionally, the B star temperature is about 0.5 the O star temperature, giving a 
number luminosity ratio of about 0.2. All in all, for photodisintegration resulting in 
PeV nucleons, the relevant photon density in the core of the Cygnus 0B2 association is 
nuv ~ 230 cm~^. The nucleus mean free path is ~ 35 kpc, corresponding to a collision 
time r = 10^ yr. Thus, the collision rate for photodisintegration in the core region is 
comparable to the hadronic interaction rate.^^ 

Since one is interested in neutrinos, it is still necessary to compare the production 
rate for charged pions in the hadronic case to the overall rate for generating neutrons. 
To assess this ratio, we made use of available high energy event simulations sho wing 
spectator nucleo n ari d pion spectra for Fe-N/p-N collisions at 10^^ and 10^^ eV |543l | 
(summarized in |544l | ) . The pion rapidity spectra in the central plateau are roughly 
energy independent, except for the widening of the plateau with energy -see |39l| . 
whereas there is a slow increase of spectator neutrons as one reaches the region of 
interest {En ^ 1 PeV). Allowing for sizeable differences in hadronic interaction m odels . 



the secondary populations are roughly 35% vr^, 45% 7, 10% nucleons, and 10% K (544 
In the energy range yielding PeV neutrons, only about 30% of the rapidity plateau 
contributes charged pions above 2 TeV. Since only half the nucleons are neutrons, we 
arrive at a ratio 

7r±(> 2 TeV) ^ 



n(~ PeV) 



(81) 



This estimate takes into account a hadronic cross section, apcp ^ ^"■'^^ w 6 x 10 cm^, and 
the generous upper limit |151j of the nucleon density ~ 30 cm~^ |153j . 
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in the hadronic interactions. 

However, photodisintegration also takes place in the outer regions of the OB associ- 
ation as long as: (i) the density of the optical photons propagating out from the core 
allows a reaction time which is smaller than the age of the cluster ~ 2.5 Myr |545| and 
(a) the diffusion front of the nuclei has passed the region in question. From Eq. ()8()|1 we 
estimate an average photon density nuv ^ 25 cm~^ out to 30 pc, which gives a reaction 
time of ~ 10^ yr. The diffusion time (~ 1.2 Myr) is a bit smaller than the age of the 
cluster, and somewhat higher than the reaction time, allowing about 90% of the nuclei 
to interact during the lifetime of the source. Thus, the production rate of neutrons 
via photodisintegration is amplified by a volume factor of 27 over the rate in the 10 
pc core. The net result of all this consideration is that the PeV neutron population is 
about an order of magnitude greater than that of the TeV charged pions 

With this in mind, we now discuss the prospects for a new multi-particle astronomy: 
neutrons as directional pointers -|- antineutrinos as inheritors of directionality. The 
basic formula that relates the neutron flux at the source {dFn/dEn) to the antineutrino 
flux observed at Earth {dFy/dEy) is ji^: 



The variables appearing in Eq. ()82|) are the antineutrino and neutron energies in the 
lab {E-jj and En)-, the antineutrino angle with respect to the direction of the neutron 
momentum, in the neutron rest-frame (^u) , and the antineutrino energy in the neutron 
rest-frame (e^). The last three variables are not observed by a laboratory neutrino- 
detector, and so are integrated over. The observable E-^j is held fixed. The delta- 
function relates the neutrino energy in the lab to the three integration variables. The 
parameters appearing in Eq. (jH^ are the neutron mass and rest-frame lifetime (m^ and 
r„), and the distance to the neutron source {D). dFn/dEn is the neutron flux at the 
source, or equivalently, the neutron flux that would be observed from the Cygnus region 
in the absence of neutron decay. Finally, dP/ de-u is the normalized probability that the 
decaying neutron in its rest-frame produces a z/g with energy e^. Setting the beta-decay 
neutrino energy ejj equal to its mean value = eo ~ (m„ — nip) [1 — m1/{mn — mp)^]/2 ^ 
0.55 MeV, we have dP/de-y = 6{e-jj — eg).^^ Here, the maximum neutrino energy in the 
neutron rest frame is Q = mn—mp—me = 0.71 MeV, and the minimum neutrino energy 



Note that = r„(ei/ + f3epCos9[^) = i?„ep(l + cos0p)/mn, where r„ — En/mn is the Lorentz 
factor, and (as usual) /3 ~ 1 is the particle's velocity in units of c. 

The delta-function in the neutron frame gives rise to a flat spectrum for the neutrino energy in the 




(82) 
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is zero in the massless limit. The expression in parentheses in Eq. ()82p is the decay 
probabihty for a neutron with lab energy En, traveling a distance D. In principle, one 
should consider a source distribution, and integrate over the volume / (PD. Instead, 
we will take D to be the 1.7 kpc distance from Earth to Cygnus 0B2; for the purpose 
of generating the associated neutrino flux, this cannot be in error by too much. 

Putting all this together, normalization to the observed "neutron" excess at ~ 
10^^ eV leads to about 20 antineutrino events at IceCube per year [i^. Note that 
anti-neutrino flux at 1 TeV may also originate in the decay of 1 PeV neutrons from 
sources whose spectrum cuts off at that energy, and hence are not subject to normal- 
ization by the anisotropy. Thus, this estimate may be regarded as very conservative. 
A direct TeV 17 ^ event in IceCube will make a showering event, which, even if seen, 
provides little angular resolution. In the energy region below 1 PeV, IceCube will 
resolve directionality only for i/^ and z/^. Fortunately, neutrino oscillations rescue the 
signal. Since the distance to the Cygnus region greatly exceeds the z/g oscillation length 
Aosc ~ 10^^(-Ep/PeV) pc (taking the solar oscillation scale 5m^ ~ 10~^eV^), the an- 
tineutrinos decohere in transit. The arriving antineutrinos are distributed over flavors, 
with the muon antineutrino flux Fp^ given by the factor ^ sin^(2^0) ~ 0.20 times the 
original Fp^ flux. The Ur flux is the same, and the Ue flux is 0.6 times the original 
flux. Here we have utilized for the solar mixing angle the most recent SNO result 
^0 ~ 32.5° [54&] . along with maximal mixing for atmospheric z/^-i^r neutrinos and a 
negligible z/g component in the third neutrino eigenstate. All in all, for a year of run- 
ning at IceCube, one conservatively expects 4 showers with energies > 1 TeV to 
cluster within 1° of the source direction, comfortably above the stated CL 

IceCube is not sensitive to TeV neutrinos from the Galactic Center, as these are above 
the IceCube horizon, where atmospheric muons will dominate over any signal. However, 
other kilometer-scale neutrino detectors, such as those planned for the Mediterranean 
Sea, may see the Galactic Center flux. 

In summary, in a few years of observation, IceCube will attain 5a sensitivity for 
discovery of the Fe— > n —>■ Vg ^ cosmic beam, providing the "smoking ice" for the 
Galactic Plane neutron hj^othesis. 



lab for fixed neutron lab-energy En — 

dP _ dcose-17 ( de-ij\ ( dP\ _ 1 

with < i;- < 2r„eo- 

The massless-neutrino approximation seems justifiable here: even an eV-mass neutrino produced 
at rest in the neutron rest-frame would have a lab energy of r„ < GeV, below threshold for 
neutrino telescopes. 
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C. Probes of new physics beyond the electroweak scale 



It is intriguing - and at the same time suggestive - that the observed flux of CRs 
bey ond the GZK-e nergy is well matched by the flux predicted for cosmogenic neutri- 



5471 IHiiH [550.]. Of course, this is not a simple coincidence: any proton flux 



nos ^ 

beyond E^^^^ is degraded in energy by photoproducing 7r° and n"^, with the latter in 
turn decaying to produce cosmogenic neutrinos. The number of neutrinos produced in 
the GZK chain reaction compensates for their lesser energy, with the result that the cos- 
mogenic flux matches well the observed CR flux beyond 10^'' eV. Recently, the prospect 
of an enhanced neutrino cross section has been explored in the context of theories with 
large compact dimensions. In these theories, the extra spatial dimensions are respon- 
sible for the extraordinary weak ness of th e gravitational force, or, in other words, the 
extreme size of the Planck mass 552 . 553^ . For example, if spacetime is taken as a di- 
rect product of a non-compact 4-dimensional manifold and a flat spatial n-torus T" (of 
common linear size 27rrc), one obtains a definite representation of this picture in which 
the effective 4-dimensional Planck scale, Mpi ~ 10^^ GeV, is related to the fundamental 
scale of gravity, M^, according to Mpj = Svr M^"^"r". Within this framework, virtual 
graviton exchange would disturb high energy neutrino interactions, and in principle, 
could increase the neutrino interaction cross section in the atmosphe re by orders o f 



magnitude beyond the SM value; namely a^N ~ [Eu/ (10^°) GeV] mb (554l . 1555 
However, it is important to stress that a cross section of ~ 100 mb would be necessary 
to obtain consistency with observed showers which start within the first 50 g/cm^ of 
the atmosphere. This is because Kaluza-Klein modes couple to neutral current s and 



the scattered neutrino carries away 90% of the incident energy per interaction |557 
Moreover, models which postulate strong neutrino interactions at super-GZK energies 
also predict that moderately penetrating showers should be produced at lower ener- 
gies, where the neutrino-nucleon cross section reaches a sub-hadronic size. Within TeV 
scale gravity a^N is likely to be sub-hadronic near the energy at which the cosmo- 
genic neut rino flux peaks, and so moderately penetrating showers should be copiously 
produced 1558^ . Certainly, the absence of moderately penetrating showers in the CR 
data sample should be understood as a serious objection to the hypothesis of neutrino 
progenitors of the super-GZK events. 

Large extra dimensions still may lead to significant increases in the neutrino cross 
section. If this scenario is true, we might hope to observe black hole (BH) produc- 
tion (somewhat more masive than Md) in elementary particle collisions with center- 



A point worth noting at this juncture: the neutrino-nucleon cross section can also be significantly 
enhanced at ce nter-of-mass energies > 100 TeV (within the SM) via electroweak instanton-induced 



processes |551l 
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of-mass energies > TeV 559l ISGOl . |56l|. In particular, BHs occurring very deep in 
the atmosphere (revealed as intermediate states of ultrahigh energy neutrino interac- 
tions) could trigger quasi - horiz ontal showers and be detected by cosmic ray observa- 
562 . 563 . 564 . 565 . 15661 - Additionally, neutrinos that traverse the atmosphere 



tones 



unscathed may produce BHs thro ugh interactions in the ice or water and be detected 
by neutrino telescopes 567, l568l |. Interestingly, cTi^n^bu oc 

fore, the non-observation of the almost guaranteed flux of cosmogenic neutrinos can 
be translated into bounds on the fundamental Planck scale. For n > 5 extra spa- 
tial dimensions compa ctified on T", recent null results from CR detectors lead to 



Md > 1.0 - 1.4 TeV [5691, Isfi. These bounds are among the most stringent and 



conservative to date. In the near futur e, PA O will provide more sensitive probes of 
TeV-scale gravity and extra dimensions [57l|. Certainly, the lack of observed deeply- 
pene trating showers can be used to place more general, model-independent, bounds on 

<yuN 



Up to now we have only discussed how to set bounds on physics beyond the SM. An 
actual discovery of new physics in cosmic rays is a tall order because of large uncer- 
tainties associated with the depth of the first interaction in the atmosphere, and the 
experimental challenges of reconstructing cosmic air showers from partial information. 
However, a similar technique to that employed in discriminating between photon and 
hadron showers can be applied to search for signatures of extra-dimensions. Specif- 
ically, if an anomalously large quasi-horizontal deep shower rate is found, it may be 
ascribed to either an enhancement of the incoming neutrino flux, or an enhancement 
in the neutrino-nucleon cross section. However, these two possibilities may be distin- 
guished by separately binning events wh ich arrive at very small angles to the horizontal, 
the so-called "Earth-skimming" events 5751 576l | . An enhanced flux will increase both 
quasi-horizontal and Earth- skimming event rates, whereas a large BH cross section 
suppresses the latter, be cause the hadronic decay products of BH evaporation do not 
escape the Earth's crust |577l | . For a more detailed discussion of neutrino interactions 
in the Earth atmosphere within TeV scale gravity scenarios see e.g. jl|. 



V. ANY GIVEN SUNDAY: COUNTDOWN TO DISCOVERY 

With data now at hand, not only there are several interesting, plausible theoretical 
models within the standard astrophysical agenda to explain all the CRs detected so far, 
but there could indeed be too many. Perhaps yet unexpected degeneracy problems will 
appear, even with the forthcoming data of the Pierre Auger Observatory, a topic which 
till now has not been a subject of debate. In our view, Occam's razor imposes that all 
standard astrophysical models be eliminated before embarking in the consequences of 
explanations involving physics beyond the standard scenarios. However, should this be 
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the case, the prospect for encountering a profound scientific revolution are endless. The 
puzzle of UHECRs may have something to say about issues as fundamental as local 
Lorentz invariance: On the one hand, the absence of photo-pion production above the 
GZK-limit would imply no cosmogenic neutrino flux and possibly undefiected pointing 
of the primary back to its source. On the other hand, a significant correlation of TeV- 
antineutrinos with directional signals at EeV energies will validate Special Relativity to 
an unprecedented boost factor of F 10^, several orders of magnitude beyond current 
limits. Additionally, contrasting the observed quasi-horizontal neutrino flux with the 
expected neutrino flux can help constrain TeV-scale gravity interactions and improve 
current bounds on the fundamental Planck scale. An optimist might even imagine the 
discovery of microscopic BHs, the telltale signature of the Universe's unseen dimensions. 
At the time of writing, new data is being collected at Pampa Amarilla. Whatever 
happens, in T-2 years we shall witness the lift-off of a new era of cosmic ray physics. 
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